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1. 


INTRODUCTION 


In  Parts  I and  II  of  this  report,  ve  have  described  our  two-, 
dimensional  computer  model  NRFL02  for  calculating  the  near-field 
external  flow  of  an  ocean  thermal  power  plant  (OTPP),  and  have  presented 
numerical  results  for  a proposed  experimental  simulation  and  for  the 
Lockheed  baseline  OTPP  design.  Here,  in  Part  III,  we  give  a listing  of 
the  computer  code  and  a brief  description  of  its  use. 

The  code  is  written  in  Fortran,  and  exists  in  both  an  IBM 
dialect  and  a CDC  dialect.  It  uses  card  input  for  five  lines  of  text, 
but  otherwise  the  input  parameters  are  specified  in  a brief  main  program 
which  must  be  re-compiled  for  every  run.  In  Sections  2 and  3 below,  we 
describe  the  input  parameters  and  the  use  of  NRFL02.  Sections  1*  and  5 
present  a brief  description  of  the  subroutines  and  their  relationships, 
and  a full  listing  of  the  code. 


r 


INPUT  PARAMETERS 


2.1  The  Main  Program 

The  parameter  values  for  NRFL02  are  not  read  as  data,  they 
are  specified  in  a brief  main  program  and  passed  to  the  subroutine  PR, 
the  master  control  program.  Thus,  in  a typical  run,  only  the  main 
program  needs  to  be  compiled.  The  main  program  is  listed  at  the 
beginning  of  Section  5 belov.  The  input  parameters  are  defined  in 
terms  of  Part  I of  this  report,  and  ve  vill  refer,  for  example,  to 
equation  (1/13)  and  page  I/lU. 

2. 2 The  Non-Dimensionalization 

The  code  NRFL02  makes  the  basic  equations  (1/13)  non- 
dimensional  by  using  the  vertical  domain  size  D (page  I/lU)  as  a 
length  scale,  scaling  the  temperature  so  that  the  ambient  temperature 


c 

Tft(z)  is  zero  at  the  bottom  of  the 

domain  and  one  at 

the  top,  and 

choosing  the  time  scale  so  that  ga 
different  variables  are,  therefore, 

is  unity.  The  units  for  the 

as  follows: 

c 

Length 

D 

(la) 

Time 

(D/goAT)3* 

(lb) 

Speed 

(DgoAT)3* 

(lc) 

Acceleration 

gaAT 

(Id) 

€ 

Diffusivity 

(D^cxAT)*5 

(le) 

Stream- Function 

( D3gaAT  J3* 

(If) 

Volume  Flux 

(D^gaAT)3* 

(lg) 

Kinetic  Energy  Density 

DgaAT 

(lh) 

where  AT  is  the  dimensional  ambient  temperature  difference  between  the 
top  and  bottom  of  the  computational  domain. 
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The  code  input  is  given  in  non-dimensional  terms,  with  the 
following  exceptions.  The  array  TFAHR(U)  of  OTPP  temperatures, 
described  in  Section  2.3,  and  the  top  and  bottom  ambient  temperatures 
TT  and  TB  (Section  2.h)  are  specified  in  degrees  Fahrenheit.  The 
vertical  domain  size  DSCALE  and  the  region  width  WIDFT  in  the  y- 
direction  (Section  2.7)  are  given  in  ft.  Note  that  for  equations  (l): 


D « 

DSCALE 

ft 

(2a) 

AT  * 

TT-TB 

ft 

(2b) 

g * 

32  ft /sec2 

• 

(2c) 

a = 

1.3  x 10~U/°T 

• 

(2d) 

These  quantities  do  not  affect  the  calculation,  but  only  the  output, 
which  is  presented  in  engineering  units,  in  terms  of  feet,  seconds, 
and  degrees  Fahrenheit. 

For  output  purposes,  z is  measured  up  from  the  top  of  the 
domain  (normally  the  water  surface).  Within  the  code,  z is  measured 
up  from  the  bottom  of  the  domain,  and  thus  z * 1 at  the  top. 

2. 3 OTPP  Design  Parameters 


The  real  arrays  AMPU,  AMPQ,  TANALF,  TFAHR 
of  dimension  U,  are  initialized  by  data  statements, 
the  dimensionless  arrays  a^,  E^,  a^,  T^,  r^,  and  z^ 
as  follows: 


, RPL,  and  ZCENT, 
and  correspond  to 
in  Section  1/3* 


s 

m 

AMPU(k) 

1 

(3a) 

=k 

m 

AMPQ(k) 

i 

(3b) 

tftnak 

m 

TANALF(k) 

• 

• 

(3c) 

Tk 

m 

(TFAHR(k) 

- TB ) / ( TT-TB ) ; 

(3d) 

rk 

9 

RPL(k) 

• 

ft 

(3e) 

9 

ZCENT(k) 

• 

ft 

(3f) 

C 
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where  z^,  being  dimensionless,  is  measured  up  from  the  domain  V>  ;tom, 
like  z. 

The  quantities  E^,  ak>  and  Tk  determine  the  inflow  values 
of  E,  w,  and  T in  equations  (l/l6),  and  are,  therefore,  irrelevant 
unless  a]£  is  positive.  Use  is  made  of  this  fact  to  provide  the  option 
of  determining  the  OTPP  outflow  temperatures  by  calculating  the  inflow 
temperatures  from  equation  (1/17)  and  adding  or  subtracting  a tempera- 
ture increment. 

The  first  two  regions  of  negative  u on  the  left  boundary, 
starting  from  the  bottom,  are  found,  and  the  corresponding  dimensionless 
mean  plant  inflow  temperatures  Tfc  are  determined  using  equation  (1/17). 
The  dimensional  plant  inflow  temperatures  TIN(l)  and  TIN(2)  are  then 
determined  using  equation  (3d),  and  TFAHR(2)  and  TFAHR ( 3 ) are  replaced 
by  the  expressions 

TFAHR(2)  * TIN ( 1 ) + TFAHR (l)  (Ua) 

TFAHR  (3 ) * TPT(2)  + TFAHR  (1*)  . (lib) 

Thus,  to  obtain  fixed  OTPP  outflow  temperatures,  positive  values  of 
al  and  a^  should  be  used,  since  TFAHR (l)  and  TFAHR (1*)  remain  unchanged. 
For  an  outflow  3°F  warmer  than  the  deepest  OTPP  inflow,  &2  should  be 
positive  and  TFAHR(l)  should  be  3;  also,  a^  should  not  be  positive, 
so  that  TFAHR(l)  is  otherwise  irrelevant.  For  an  outflow  3°F  cooler 
than  the  second  deepest  OTPP  inflow,  a^  should  be  positive,  TFAHR (1*) 
should  be  -3,  and  a^  should  not  be  positive.  If  there  is  only  one 
OTPP  inflow  region,  then  TIN(2)  will  be  zero,  and  the  original 
TFAHR(3)  value  will  be  replaced  by  TFAHR(U). 

It  has  been  suggested  that  the  warm  and  cold  outflows  from 
the  plant  should  be  mixed.  This  can  be  accomplished  consistently  by 
“•king  *k  Posit^ve  for  any  two  k values,  and  making  the  corresponding 
z^  and  r^  values  equal. 
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The  Ambient  Ocean 


2.k 


The  dimensionless  ambient  temperature  distribution  is  given 
by  equation  (l/l8a),  where 


-ZTC 

(5a) 

WID 

(5b) 

and  T and  T are  defined  so  that  T (0)  = 0 and  T (l)  = 1.  The 
**  t a a 

dimensionless  ambient  turbulent  kinetic  energy  distribution  is  given 
by 

E&(z)  = QAMB  exp  | (z  - 1)/QDEP j . (6) 

The  dimensionless  parameters  QAMB  and  QDEP  correspond  to  E and  z_ 

o E 

in  equation  (l/l8c). 

2. 5 The  Turbulence  Model 


The  parameters  defining  the  turbulence  model  are: 


L = 

TLEN 

9 

(7a) 

Cf  " 

TCOF 

9 

(7b) 

cs  * 

CSTRAT 

9 

(7c) 

Cw  " 

CW 

• 

(7d) 

Here  L and  TLEN  are,  of  course,  dimensionless;  the  dimensional 
turbulence  length  scale  is  TLEN  x DSCALE  ft. 

In  addition,  the  parameter  PORC  is  set  in  this  section  of 
the  main  program.  PORC  replaces  the  constant  ^ in  equation  (l/22b), 
determining  the  porosity  distribution  which  stops  the  reflection  of 
internal  waves  by  the  right-hand  boundary. 


2. 6 The  Computational  Mesh  and  Time-Step 

The  number  of  mesh  points  in  the  vertical  direction  is  input 
as  the  parameter  JJ  in  this  section;  JJ  must  he  no  greater  than  the 
vertical  dimension  in  the  array  declarations  in  subroutines  PR  and 
MARCH  (1*0  in  the  listing  at  card  B28  in  Section  5,  but  see  Section  3.2), 
The  number  of  mesh  points  in  the  horizontal  direction  is  II,  and  is 
set  in  the  subroutine  PR. 

The  computational  domain  is  rectangular,  with 

0 < z < 1 , (8a) 

0 5 x £ y , (8b) 

in  dimensionless  terms,  cf.  page  l/lh.  Here 

Y = WIDTH  . (9) 

The  mesh  intervals  6x  and  6z  increase  in  geometric  progression,  with 
the  rate  of  increase  defined  by  the  input  parameters 


XSTRCH  = log(6xright/fixleft) 


ZSTRCH  = log(6z,  ..  / 6z , ) 

bottom  top 


The  total  dimensionless  calculation  time  is 


(10a) 


(10b) 


t , = TOTIME  x WIDTH/u 

end  max 


where  u is  the  maximum  of  the  imposed  horizontal  velocity  on  the 
max 

left  boundary.  The  time-step  is  determined  from  equation  (1/3^)*  vith 
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2.7 


The  Output  Options 


The  code  prints  a diagnostic  line  of  maxima  and  integrals, 
for  every  time  step.  In  addition,  it  prints  a variety  of  contour 
plots  and  curves  NOUTPT  times  per  run,  including  the  final  print  which 
gives  greater  detail.  Afterwards,  NCOPYS  copies  of  a result  summary 
are  printed;  this  result  summary  is  illustrated  by  Figures  1/3  and  I A. 

DSCALE  is  the  depth  D of  the  computational  region,  in  feet. 
WIDFT  is  the  width  of  the  region,  in  the  y-direction.  It  is  multiplied 
by  the  dimensional  stream-function  to  give  the  volume  flux  in  ft^/sec. 


3. 


USE  OF  NRFL02 


3*1  Choice  of  Input  Parameters 

Once  a particular  dimensional  temperature  profile  (l/l8a) 
and  domain  size  D = DSCALE  have  been  chosen,  then  TT  and  TB  can  be 
calculated,  and  the  units  (l)  can  be  found.  Because  the  stratification 
tends  to  produce  horizontal  flow,  large  D values  are  not  required. 

Since  the  external  flow  for  the  OTPP  designs  at  present  under 
consideration  is  not  two-dimensional,  the  determination  of  the  design 
parameters  (Section  2.3)  to  be  used  in  the  model  is  difficult.  An 
example  is  given  in  Part  II  of  this  report  for  the  Lockheed  baseline 
design;  the  greatest  uncertainty  is  in  the  appropriate  value  of 
WIDFT.  The  arrays  AKPU,  RPL,  and  ZCENT  must,  of  course,  be  scaled. 

The  ambient  turbulence  must  be  small  enough  to  be  innocuous, 

-3  —6 

but  must  not  be  zero.  Dimensionless  QAMB  values  between  10~  and  10~ 
have  been  used,  with  QDEP  from  ~ to  10.  The  results  are  not  at  all 
sensitive  to  these  values. 

The  turbulence  model  parameters  are  a major  problem.  Our 
present  choices 


TLEN 

= L * 

RPL  x 0.1+  , 

(13a) 

TCOF 

II 

V. 

o 

II 

0.5 

(13b) 

CSTRAT 

II 

</) 

o 

II 

0.1 

(13c) 

CW 

Cv 

0.1  , 

(13d) 

are  based  on  our  experience  in  the  simulation  of  turbulent  submarine 
wakes.  These  parameters,  and  the  corresponding  parameters  in  NRFL02 
successors,  must  be  tuned  to  give  the  best  agreement  with  a wide  range 
of  experimental  observations.  The  porosity  coefficient  TCOF  should  be 
about  it  is  a part  of  the  numerical  method,  not  the  turbulence  model. 


The  domain  WIDTH  (9)  should  be  between  2 and  U.  The  calcu- 
lations in  Part  I,  with  the  maximum  possible  port  separation,  used 
Y - 2 for  display  convenience;  a slightly  larger  value  would  have  given 
a better  far-field  result.  Values  of  XSTRCH  from  0.5  to  2 are  appro- 
priate; the  calculation  needs  good  horizontal  resolution  only  near  the 
left  boundary.  The  number  II  of  mesh  points  in  the  x-direction  (see 
Section  3.2)  should  be  such  that  the  mesh  interval  6x  at  the  left  is 
about  half  of  the  smallest  RPL  value.  The  quantities  JJ  and  ZSTRCH 
should  similarly  be  chosen  so  that  the  mesh  interval  6z  is  about  one 
third  of  any  neighboring  RPL  values.  Less  resolution  is  needed  in  z 
ranges  well  away  from  any  significant  ZCENT  values,  and  ZSTRCH  magni- 
tudes of  1 or  more  can  sometimes  be  specified. 

TOTIME  values  of  10  to  20  have  been  used  in  our  computations. 
If  the  solutions  do  not  become  effectively  steady  in  this  time,  then 
the  oscillations  appear  to  continue  indefinitely;  our  problems  with 
instability  of  the  mean  flow  solutions  have  not  yet  been  entirely 
solved.  The  time-step  parameter  CFL  should  be  about  0.6;  larger 
values  occasionally  result  in  numerical  instability. 

3- 2 Dimension  Statements 


The  computational  arrays  are  dimensioned  in  the  subroutines 
PP  and  MARCH,  in  the  statements  labeled  B2l*  to  B^3  and  D22  to  D37 
in  the  listing  in  Section  5*  The  value  of  JJ,  specified  in  the  main 
program,  can  be  any  number  up  to  the  corresponding  dimension  in  the 
arrays  [ho  in  the  listing).  The  value  of  II  must  agree  with  the  first 
dimension  of  the  two-dimensional  array  declaration  (12  in  statements 
B28  to  B31  and  D22  to  D25);  therefore,  II  is  set  to  12  in  statement 
B60.  Values  of  II  up  to  25  can  be  used  by  changing  the  declarations 
of  the  two-dimensional  arrays;  larger  values  require  changing  the 
dimension  statements  for  the  one-dimensional  arrays  associated  with 
the  x-direction. 


3.3 


Input  and  Temporary  Storage  Streams 


The  required  streams  are  listed  in  CDC  form  in  the  comment 
card  at  the  beginning  of  the  main  program.  Five  cards  of  text  (up  to 
60  characters  per  card)  are  read  from  stream  5 and  printed  on  the 
output;  the  first  card  is  printed  as  a heading  for  all  the  plots. 

The  streams  8 to  13  are  used  for  temporary  storage,  with  a 
printer  stream  format.  Some  of  the  output  routines  write  to  these 
streams,  and  the  subroutine  COMBIN  then  reads  from  specified  streams, 
combines  several  streams  as  needed,  and  copies  them  to  the  printer 
stream. 

The  printer  output  is  on  stream  6.  A printer  line  of  up  to 
128  characters,  plus  the  control  character,  is  produced. 

3- U Output 

The  text  from  the  five  input  cards  is  printed  first,  followed 
by  a list  of  the  input  parameter  values  from  the  main  program.  This 
section  of  text  is  repeated  in  the  result  summary  output , as  at  the  bottom 
of  Figures  1/3  to  1/8.  The  non-uniform  meshes  and  intervals,  and  the 
tine-step,  the  number  of  time-steps  and  the  total  dimensionless  time, 
determined  by  the  code,  are  also  printed. 

For  every  pair  of  leapfrog  time-steps,  the  code  outputs  a 
diagnostic  printer  line  giving  various  maxima  and  minima,  together  with 
the  Fahrenheit  OTPP  inflow  temperatures  determined  using  equation  (1/17) 
and  equation  (3d).  This  output  can  be  used  to  investigate  convergence 
to  a steady  solution,  as  well  as  for  detecting  errors  or  instability. 
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The  total  number  of  time-steps  is  split  into  NOUTPT  equal 
parts,  and  after  each  section  a number  of  printer  plots  ewe  generated. 
The  first  three  show  the  volume  flux  added  to  the  ocean  below  depth  z 
on  the  left  boundary  and  on  the  right  boundary,  and  below  temperature  T 
on  the  right  boundary.  This  output  is  illustrated,  for  the  strong- flow 
case  of  Part  I,  in  Figure  1. 

A plot  of  the  surface  temperature,  as  a function  of  x,  is 
presented  next,  showing  the  temperature  loss  through  turbulent  mixing 
as  the  surface  layers  ewe  drawn  towards  the  OTPP  inflow. 

Contour  plots  are  printed  next,  displaying  seven  vewiables 
as  functions  of  x and  z.  The  variables  are  as  follows: 


Turbulent  kinetic  energy  density 


E(ft2/sec2) 

J 

(lUa) 

Turbulent  diffusivity  K(ft2/sec) 

* 

(lUb) 

Diffusivity  factor  cgN2L2/E 

> 

(Ike) 

2 i — 2 

Richardson  number  N /-^u,  + u.  . ) 

* 9 J J * 

J 

(Ikd) 

Reynolds  number 

{2*Ui,J  + uJ,i)  ^/K 

» 

(lUe) 

Temperature  T(°F) 

» 

(l^f) 

Volume  flux  function 

WIDFT  x ^(ft3/sec) 

• 

(lUg) 

These  ewe  contoured  over  the  whole  domain,  with  a specification  of  12 
contour  intervals.  At  the  end  of  the  run,  the  seven  plots  are  repeated 
for  Just  the  left  half  of  the  domain,  double  scale,  using  20  contour 
intervals.  This  gives  higher  resolution  in  the  region  where  smaller 
length  scales  are  present. 


lfc 


FIGURE  1.  Illustration  of  the  Environmental  Impact  Plots 


The  three  printer  plots  show  the  volume  flux  in 
ft  /sec  added  to  the  "ocean"  below  depth  z ft  on 
the  left  and  right  boundaries,  and  below  tempera- 
ture T on  the  right  boundary,  for  the  NRFL02  simu- 
lation of  the  proposed  strong-flow  experiment 
reported  in  Part  I.  On  the  left  boundary,  a flux 
1.U1*  ft3/sec  (for  1 ft  of  tank  width)  is  added  in 
the  bottom  1 ft , and  an  equal  flux  is  removed  in 
the  1 ft  interval  at  the  surface.  As  a result  of 
the  turbulent  transport,  entrainment,  and  recir- 
culation in  the  near-field,  the  far- field  impact 
is  the  addition  of  0.606  ft3/sec  below  a depth  of 
1.9  feet  in  the  temperature  range  from  1*3°F  to 
59°F , and  the  removal  of  an  equal  flux  from  the 
shallower  water,  between  59°F  and  80°F. 


’ 

Finally , the  result*  are  summarized  using  the  format  of 
Figures  1/3  and  I A.  One  copy  of  this  format  is  generated  at  the 
intermediate  stages,  and  NCOPYS  copies  are  generated  at  the  end,  so 
that  a copy  can  be  photographed,  and  so  that  copies  can  be  sent  to 
OTPP  designers  direct  from  the  printer. 
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1*.  BRIEF  DESCRIPTION  OF  THE  SUBROUTINES 

fc.l  NRFL02  and  PR 


In  the  program  listing  in  Section  5,  each  subroutine  baa  a 
different  alphabet  identifier.  In  this  section,  we  give  the  identifier 
in  brackets  (note  that  this  identifier  is  not  an  argument).  Vhere 
appropriate,  we  refer  also  to  the  individual  card  identifier. 

The  main  program  NRFL02(A)  must  be  compiled  for  every  run. 

It  calls  PR(B) , the  master  control  program,  which  does  all  the  work. 

PR  and  the  subroutine  MARCH(D)  must  be  recompiled  every  time  II,  the 
number  of  mesh  points  in  the  x-direction,  is  changed  (see  Section  3.2). 

The  other  subroutines  need  only  be  compiled  once. 

The  master  control  program  PR  receives  the  input  parameters 
from  the  main  program  (some  under  different  names,  because  they  are 
passed  elsewhere  in  common  blocks).  It  reads  the  text  cards,  finds 
the  date,  and  writes  (B75)  the  text  and  the  list  of  input  parameter 
values  to  the  printer  and  to  stream  10  (for  later  printing).  It  sets 
up  the  mesh  (3100)  and  the  boundary  and  initial  conditions  (B158),  with 
the  subprogram  FUNC(C)  used  for  the  function  f(x)  in  equation  ( l/lh ). 

It  calculates  the  time-step  and  the  time-step  counts  (B230),  and  sets 
the  auxiliary  arrays  for  the  calculation  (B2U8)  and  the  arrays  for 
solving  Poisson's  equation  ( 1/31 ) [using  the  subroutine  ROPT(l)]. 

The  leapfrog  time-stepping  is  done  by  the  subroutine  MARCH(D). 
There  are  two  arrays  for  each  of  the  four  variables  u,  w,  T,  and  E, 
one  for  the  even  time-steps  and  one  for  the  odd.  MARCH  is  called  twice 
for  each  pair  of  time-steps,  to  replace  first  one  set  of  arrays  and  then 
the  other  set  with  its  new  values.  The  main  loop  extends  from  B293  to 
the  end,  with  the  initialization  and  parasite  removal  (page  1/27) 
followed  by  the  double  call  to  MARCH,  the  call  to  OUTPUT(M)  which  controls 
all  the  printing,  and  the  modification  of  the  plant  outflow  temperatures 
according  to  equations  (U)  and  (3d). 
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^*•2  MARCH  and  Associated  Subroutines 

The  subroutine  MARCH(D)  extrapolates  the  pressure  for 
equation  (1/29)  and  evaluates  the  turbulent  diffusivity  (Dl*5).  It 
prepares  arrays  classifying  the  mesh  points  on  the  left  and  right 
boundaries  as  inflow,  slow  outflow,  or  fast  outflow  (D6l),  for  applying 
the  boundary  conditions  in  Sections  1/3  and  1/1*  and  the  prescription 
on  page  I/2l*  for  fast  outflow. 

The  subroutines  FCSTD(E),  FCSTU(F),  and  FSCTW(G)  advance 
E,  T,  u,  and  w according  to  equations  (1/26),  see  card  D102.  These 
subroutines  use  TRIDl(H)  for  the  tridiagonal  implicit  solutions. 

The  pressure  operations  described  on  page  (1/26)  are  applied  starting 
with  card  Dill*,  using  the  subroutines  TSTEP(K)  for  equations  (1/30) 
and  POISB(J)  for  the  Poisson  equation  (1/31).  The  subroutine  BOUND(L) 
is  then  called,  to  apply  the  boundary  conditions  on  the  top  and  sides 
using  the  arrays  determined  earlier.  This  completes  the  time-step. 

^ • 3 OUTPUT  and  Associated  Subroutines 

The  subroutine  OUTPUT(M)  is  called  from  PR  every  time-step 
pair,  and  generates  the  output  described  in  Section  3.1*.  It  evaluates 
the  various  maxima  and  minima  (M29),  the  divergence  (M51),  and  the 
inflow  temperatures  (M59),  and  prints  the  diagnostic  line.  Control  is 
then  returned  to  PR  unless  the  count  is  satisfied  for  further  output 
(M93). 

The  inflow  temperatures  are  written  to  the  end  of  stream  10, 
for  later  printing  with  the  results  summary  (M100).  The  mesh  arrays 
are  made  dimensional  (Mill*)  and  the  speed  scale  (lc)  is  found  (M130). 

The  Figure  1 environmental  impact  plots  and  the  surface  temperature 
plot  are  printed  (M129),  using  the  subroutine  XPLOT(R),  which  calls 
MAXMIN(S)  and  SH0RT7(T).  The  seven  contour  plots  are  printed  in  the 
section  from  M150  to  M219,  using  the  subroutine  CONTOR(N).  The  loop 
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starting  at  M153  produces  two  series  of  seven  plots  at  the  end  of  the 
calculation.  The  second  series  displays  only  the  left-half  of  the 
computational  region,  with  more  contours  and  with  twice  the  resolution 
of  the  first  series.  The  subroutine  BD(P)  is  used  to  prepare  some  of 
the  seven  functions  for  CONTOR,  by  applying  symmetry  boundary  conditions 
to  specify  undetermined  boundary  values. 

The  plots  in  the  results  summary  of  T and  u at  the  left  and 
right  boundaries,  together  with  the  ambient  temperature,  are  produced 
by  the  subroutine  VPIX)T(Q).  The  speed  and  temperature  plotting  ranges 
are  evaluated  after  card  M222,  using  the  subroutine  MAXMIN(S).  All 
four  plots  in  the  result  summary,  as  shown  in  Figures  1/3  and  iA,  use 
the  same  call  M263,  with  a loop  M2U5  for  left  and  right  and  a loop  M259 
for  small  and  large.  The  plots  are  written  on  streams  8,  9,  11,  and  12; 
the  scales  are  then  written  beneath  the  plots  (M26M. 

The  results  summary  is  written  starting  at  M279;  if  multiple 
copies  are  appropriate,  it  is  written  to  stream  13  (M28l)  and  NCOPYS 
copies  are  written  later  to  stream  6 (M3^). 

First,  the  temperature  and  the  volume  flux  function  (lUg) 
are  contoured  again,  this  time  with  a specification  of  9 contour 
intervals  (M285).  Then,  the  small  size  VPLOT  output  on  streams  8 and  9 
is  combined  with  the  results  summary  on  stream  10,  using  the  subroutine 
COMBIN(O),  to  generate  the  format  in  the  lower  panel  of  Figure  l/3 
(M308).  Next,  the  diffusivity  is  contoured  (M313).  The  large-scale 
VPLOT  output  on  streams  11  and  12  is  then  combined,  as  in  the  center 
panel  of  Figure  I A (M323).  Finally,  the  results  summary  on  stream 
10  is  split  up  between  streams  8 and  9 (M325)  and  recombined  (M3^2) 
to  produce  the  format  in  the  lower  panel  of  Figure  I A. 

This  completes  the  results  summary,  and  after  restoring 
the  dimensionless  meshes  (M353)  control  is  returned  to  PR. 
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5.  , NRFL02  LISTING 

f 

The  code  listing  is  reproduced  photographically  on  the 
following  pages,  in  an  IBM  dialect.  Each  card  has  an  identifier, 
the  letter  being  different  for  each  subroutine.  The  subroutines 
C are  listed  below,  with  their  identifiers,  names  of  the  calling 

subroutines,  and  the  purpose. 


( 

IDENTIFIER 

AND  NAME 

CALLED  FROM 

PURPOSE  ■ 

A 

NRFL02 

Initialize  parameters  for  PR 

B 

PR 

NRFL02 

Master  control  program 

C 

FUNC 

PR 

Function  f(x)  in  equation  ( 1/lU ) 

D 

MARCH 

PR 

Advance  one  leapfrog  time-step 

E 

FCSTD 

MARCH 

Advance  E and  T 

F 

FCSTU 

MARCH 

Advance  u 

f 

G 

FCSTW 

MARCH 

Advance  w 

H 

TRID1 

FCSTD 

FCSTU 

Tridiagonal  implicit  solution 

FCSTW 

c 

I 

ROPT 

PR 

Set  arrays  for  POISB 

J 

POISB 

MARCH 

Poisson  solver 

K 

TSTEP 

MARCH 

Add  pressure  gradient  to  flow 

L 

BOUND 

MARCH 

Apply  boundary  conditions  i 

c 

M 

OUTPUT 

PR 

Output  control  program 

N 

CONTOR 

OUTPUT 

Contour  program 

0 

COMBIN 

OUTPUT 

Combine  output  streams 

P 

BD 

OUTPUT 

Prepare  array  boundaries  for 

c 

CONTOR 

Q 

VP  LOT 

OUTPUT 

Vertical  function  plots 

• 

R 

XPLOT 

OUTPUT 

Horizontal  function  plots 

S 

MAXMIN 

OUTPUT 

Maximum  and  minimum 

XPLOT 

•4 

1 

T 

SH0RT7 

XPLOT 

Short  format  for  numbers 

i * 

L 
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PROGRAM  NRFL02 1 1 NPUT. OUTPUT, TAPES" INPUT, TAPE 6"0UT PUT , TAPES, TAPER, 
I TAPEiO) 

INITIALIZE  ALL  PARAMETERS  ANO  CONSTANTS  ••• 

IMPLICIT  REAL'S! A-H.O-Z) 

• • • • OCEAN  ENGINE  PARAMETERS  • • * 

REAL  A IP  J(  4)  , AMPQ!  A I , T A NAL  F(  4)  ,TF  AMR  ( 4)  , RPL  ( 4 1 , ICE  MT 1 4 I 
DATA  AM»U/-.12,.12,2*0./ 

UATA  AAPQ/O. ,.2,2*0./ 

DATA  T ANALF/4*0. / 

DATA  T*4MR/-3.,70.,2*0./ 

OATA  R»L/4*.2/ 

OATA  ZC5NT/.3, .4,2*0./ 


AMBIENT  OCEAN 
TT"80 
TB-6T 
ZTC— .? 
MID".  I 
OAMS-.OOl 
MOEP-.l 


TURBULENCE  MOOEL  PARAMETERS 
TLEN". 03 
TCOF-.l 
CSTRAT..1 
CM". I 
PURC". i 


ME  SH 

XSTACM-LOSI RIGHT  DX  / LEFT  DX  I 
ZSTRCM-LOGI  iOTTON  01  / TOP  OZ  » 

JJ-12 

W10TH-2 

xstr;h«i 

ZSTRCMO 


TIME  AMI)  STE* 

TOT  I ME  • 1 0 
CFL-.6 


OUTPUT 

NOUTPT  "2 
DSC  AL  5" ISO 
MIDFT "l JO 
NCO*TS«i 

call  pr  i jj, ampu,ampo,tanalf,tfahr,apl,z:ent,tt,tb,ztc,mio,3amb, 
iqdep,tle«.tcof,cstrat,cm,poac,ni'otm,xstr:h, zstrch.totime.cf, ,noutp 
2 T , use  AL  E , A I DF  T , .NCOPr  S I 
STOP 
END 


SUBROUTINE  PR  UJT,AMPJ,AMPO,TANALF,TFAHR,RPL,ZCENT,TT,TB.ZTC,MI  B 
lO.OAMB, J OEP, TL EN, TC3F, TC S T, T CM, PORC, Ml OTH,XSTRCM,ZSTRCH,TOTI ME,  CFL  B 
2 , NJUT  PT  , 3>CALE«m  IOFT  .NCOPTSI  B 

•••  MASTER  CONTROL  PROGRAM  •• • B 

IMPLICIT  real*bi a-n.o-zi  b 

common  /ni/  ot.time  b 

COMMON  /N2/  C0F,RLAM,P1  , GA  P,  DE P,  CSTR AT  ,CM  B 

COMMON  /N3/  II,JJ,tU,JJl,  II2.JJ2  B 

COMMON  / N4/  AAR, AAZ.XAPPA  B 

COMMON  /NT/  A1 , A2, A3, Bl, M2, MS, MI ,MK  B 

COMMON  /SI  2/  F,E,A0V,AV2,S,P,*E  B 

COMMON  /Nl 3/  ALP,BET,GA",DEN,N,SO,ZS,S01,S02,ZSI,ZS2,ZFTU,ZFTM,Z  B 
1 F 0 J ,ZFOA,XFTU,  XFTM.XFUU.XFOM.XPOR  B 

COMMON  /N 14/  0X1 ,0X2 , OXP ,0 XM , OXT , OXG ,0X0 , 0XU, DXM ,0Z 1 , DZ2 , 0 ZP , 0 ZM  B 
1,OZT,OZG,OZO,OZU.OZM,ZBL,Z3R,ZML.ZWR  B 

COMMON  /N1S/  A ,X,Z , V B 

COMMON  /N16/  B ,BF, UBNO, F AC  B 

REAL  ASS,  SORT  ■ B 

REAL  T,ZT,Q,ZQ,F,E, A0V, A V2  B 
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REAL  vJ,ZU,W,ZM,S,P,PE 

B 

22 

c 

B 

23 

c 

OIMENSIOI  STATEMENTS  IN  MAIN  A NO 

MARCH 

B 

24 

c 

10  ARRAYS  CA i BE  LEFT  WITH  LARGER  OIMENSION,  BUT  FIRST  DIMENSION 

B 

25 

c 

OF  20  array;  MUST  BE  RIGHT. 

B 

26 

c 

B 

27 

DIMENSION  UI12.40),  2UI  12.40). 

W 1 12 , 40 1 , 2 Ml  12 , 40) 

, SI  12,40) , 

PI 

B 

28 

112.40).  Pc  I 12. 40 ) 

B 

29 

DIMENSION  T I 12.40) • 2TI 12,40). 

0(12,40),  20(12,40) 

, F(  12,40) , 

El 

B 

30 

112,40),  W2I12.40),  ADVI 1 2 ,40) 

B 

31 

DIMENSION  RI2S).  XI2S),  2140), 

YI 40) 

B 

32 

DIMENSION  OR  1(25),  0X2(23),  0XPI2S),  0XMI25I,  DXTI25),  0X3(25) 

t 

B 

33 

10X0(25),  0XUI25),  0XWI25) 

B 

34 

DIMENSION  021(40).  022(40),  02PI40),  02M(40),  02T(40),  DZ3I40) 

9 

B 

35 

1020(401,  02UI43I,  D2WI4D),  23LI40),  ZBR(40),  2WLI40) 

, 2MR ( 43 ) 

6 

36 

DIMENSION  A1I25),  A2I25),  A 3 ( 2 S ) , Blt40),  B2I40), 

B3 ( 40 ) , W I ( 25 ) 

B 

37 

1,  MKI40) 

B 

3B 

DIMENSION  B(  40 ) , FACI40),  BFI40) 

B 

39 

DIMENSION  UBN0I40) 

B 

40 

DIMENSION  ALPI40I.  BETI40),  GAM(40),  OENI40),  H(40),  S0I4D),  Z SI 

B 

41 

140),  SOU  40)  , S02  ( 43  ) , 2S1I40), 

2S2I40I,  XFTUI40),  XFTM(AO),  XFOUI 

B 

42 

240),  XFUM40),  2FTUI40),  2FTMI43),  2F0UI40),  ZFQM(43),  XP0RI40) 

B 

43 

c 

B 

44 

c 

B 

45 

REAL  A<PJ(4),AMPQ(4),TAN41F(4I 

.TFAHRI4) ,RPL(4),2CENT(4) 

B 

46 

COMMON  /TXTC/  ITxri60),n»Y.JTXT(60) 

B 

47 

DIMENSION  TINI2) 

B 

48 

c 

B 

49 

c 

8 

50 

c 

TXT  IS  TEXT  FOR  EACH  RUN,  ON  XPLOT  AND  CONTOR  ANO  VPLOT 

OUTPUT. 

B 

51 

c 

READ  OFF  DATA  CARO,  (60A1)  . 

B 

52 

RE AO  (5,10)  I TXT , JTXT 

B 

53 

10 

FORMAT  ( 60 A1 / ( 15A4I) 

B 

54 

CALL  I 3 AY  (DAY) 

B 

55 

c 

B 

56 

11-12 

B 

57 

JJ-20 

B 

58 

I 1-8 

B 

59 

11-12 

B 

60 

JJ-JJT 

B 

61 

TIME-TOTIME 

R 

62 

GAP -WIDTH 

B 

63 

rlam-tlen 

B 

64 

COF-TCOF 

B 

65 

CSTRAT-TCST 

B 

66 

CM-TCM 

B 

67 

c 

B 

6 B 

c 

ENO  INITIALIZATIONS 

e 

69 

c 

B 

70 

c 

B 

71 

c 

* * * 

B 

72 

c 

B 

73 

00  20  1-6,10,4 

B 

74 

20 

WRITE  (1.33)  I TXT, JTXT, AMP U, AMPQ, T ANAL F , TF AHR , R PL . 

2CENT.TT ,TB, 

zr 

B 

75 

IC,MID,0ARB,Q0EP,RLAM,C0F,CST RAT, 

CM, P ORC, GAP, XSTRCM.ZSIRCH, II ,JJ, 

T I 

B 

76 

2ME,CFL,NOUTPT,OSCALE,WIOFT,UAY 

B 

77 

JO 

FORMAT  (1H1//1X60A1//4I  1R,  15 A4 /)/////•  AMPU  *. 

4F 12.4/'  AMPQ 

B 

78 

1 * , 4F  1 2 . 3/  ' TANALF  •.4F12.3/' 

TF AHR  *,4F12.l/'  RPL  ',4F12 

.3 

6 

79 

2/*  2CENT  '.4F12.3/*  TT,TB,2TC, 

MIO  * ,F4.1,F  12.  1,F  14 

. 3.F12.J,/* 

OA 

B 

80 

JMB.03EP  *,1P,024.1,3P,F9.2/*  TLE  N,  CF  ,C  S , C W , » ORC  • . 5F7 

.3,/*  GAP,  XSTR 

B 

81 

4CH,ZSTRCH,II,JJ* ,3F8. 3,214/*  T 1 ME ,CF L .NOUTPT , DSC ALE. 

MIOFT* ,F6.l, 

F & 

B 

82 

5.2, 13, 2FT.1/ IX, AS, ////////) 

B 

83 

r 

B 

84 

c 

R 

65 

CALL  INOUMP 

B 

86 

CALL  N3UN0F 

B 

87 

PI-4-DATANI1 .000 ) 

B 

88 

KAPPA-2 

B 

89 

KAPPA- 3 

B 

90 

IIl«tl-l 

B 

91 

112-11*2 

B 

92 

JJ1-JJ-1 

B 

93 

JJ2-JJ-2 

B 

94 

1 2- 1 1 2 

B 

95 

J1-JJ1 

B 

96 

c 

B 

97 

c 

B 

98 

c 

SET  UP  MESH 

B 

99 

FMJLT-l-XSTRCH/I 12 

B 

100 

2MJLT-1-2STRCH/JJ2 

6 

101 

23 


OZ-1 

B 

102 

OX-DZ 

8 

103 

DEP-1 

6 

106 

c 

B 

105 

c 

B 

106 

RFM-FNJLT44.5D0 

B 

107 

oxiiii*dx/fmult 

B 

108 

0X211  l-OX/RFN 

8 

109 

xm-o 

B 

110 

Rill  — OX/RFM/I  14RFN) 

B 

111 

00  40  1-2,11 

B 

112 

0X2(1 l*QX2(I“ll-FMULT 

B 

113 

0X11  N-OXKl-1  )*FMULT 

B 

116 

xi  n«xii-i)»oxi<  n 

B 

115 

43 

R( I I-R ( 1-1 )*0X2< I-ll 

B 

116 

FACT-GAP/XUUI 

B 

117 

00  53  I-l.II 

B 

118 

X( I )*X< I l-FACT 

B 

119 

R<I!»R< 1 l-FACT 

B 

120 

0X1  ( 1 I«DX1 ( I l-FACT 

B 

121 

53 

0X2(  I » *0X2  ( 1 l-FACT 

B 

122 

PRINT  40 

B 

123 

63 

F OR H AT  (//•  OX  1,0X2, X,R*  ,/ I 

B 

126 

PRINT  TO,  <0X1111,1-1, III 

B 

125 

PRINT  TO,  (0X21 I>, 1-1,11 > 

B 

126 

PRINT  TO,  IXII),I>1,III 

B 

127 

PRINT  TO,  (RUI,I-1,UI 

B 

128 

TO 

FORMAT  1 / * 10  FI  2.  5)) 

B 

129 

C 

B 

130 

c 

B 

131 

RFWULT--.5D0 

B 

132 

DZim  -02/zhult 

B 

133 

0221  1 1 » 3Z/RFM 

B 

136 

Y<1  1*0 

B 

135 

1(11- 32/RFM/l  1-RFM) 

B 

136 

30  80  J*  2,  OJ 

B 

137 

0 2 1 ( J)*DZ1( J-l l-ZMULT 

B 

138 

3Z2<  J) *0Z2( J-l 1-ZHULT 

B 

139 

V(J»-YCJ-i»*OZ 1C  J» 

B 

160 

Z< J)«Z< J-1I-0Z2I J-ll 

B 

161 

80 

CONTINUE 

B 

162 

FACT* 1/Y( J J 1 1 

B 

163 

00  93  J-l, 00 

B 

166 

3Z1!  J)«DZl( J l-FACT 

B 

145 

0Z2<  JI«0Z2<  J l-FACT 

B 

146 

T(  J)*V( JI-FACT 

B 

147 

90 

Z(J)*Z( JI-FACT 

B 

146 

c 

B 

149 

PRINT  100 

B 

150 

130 

FORMAT  <//•  0Zl,0Z2»Y,Z',/l 

B 

151 

PRINT  TO,  (0Z1  ( J ) , J-l , J J t 

B 

152 

PRINT  TO,  I0Z2I 0>,J*1,J0) 

B 

153 

PRINT  T3,  < T ( J ) ,J*1,JJI 

B 

154 

PRINT  TO,  (Z(OI,J*l,JOI 

B 

155 

c 

B 

156 

c 

6 

1ST 

z • 

0 0 0 B 

158 

c 

B 

159 

C LEFT 

AN9  INITIAL  U IS  UBN3IJI 

B 

160 

C LEFT 

U,  0,  T,  PERMANENT  IN  ARRAYS,  WHERE  U.GT.O 

B 

161 

C RIGHT 

AND  INITIAL  T IS  B <J> 

B 

162 

C RIGHT 

AN3  INITIAL  0 IS  BFIJI 

B 

163 

C RIGHT 

AN3  INITIAL  W IS  ZERO 

B 

164 

C 

B 

165 

0AR5T * I ATANI  < 0 EP-Z TC 1/ HI  01 

B 

166 

DARGB*  3AT  ANI  ( 0 ZTC  I / HI  0 1 

B 

167 

ZERO-O 

6 

168 

UMAX— 0 

B 

169 

03  120  J-l.JJ 

B 

170 

OJNDI JI-0 

B 

171 

W(1,JI>3 

B 

172 

4-1.3-10 

B 

173 

C-0 

8 

174 

0-0 

B 

175 

33  110  K-1,A 

B 

176 

AA«FUNC(IZ(J»-ZCENT(XI|/RPL(XI  1 

B 

177 

AB-FUNCI  (Y<  JI-ZCENTIK)  1 / RPlUI  1 

B 

178 

A'J*  AMP J ( K I /R  PL  (X» 

B 

179 

AV*044  XII AU, ZERO) 

B 

180 

U8N0C  J I-U3N0I J I ♦ AJ  • AA 

B 

181 

f 


f 


c 


c 


c 


c 


c 


c 


c 


l 


A-A»AV*AA  B 

C»C44V»AV*AV-AMPQ{  M-AA  B 

0-3*4V-(TFAHR( M-TBI/ITT-T  8)  *44  B 

WCl,J|-WClfJI*AV*AB*TANALF CM  B 

110  CONTINUE  B 

QC1.J1-C/A  B 

TCl,J)-0/A  B 

Ul 1. JI-UBNDI J)  B 

UNAX- ON AX 1 ( UMAX,  OABSI  U3N0(  J)  ) I B 

8IJI-I0ATANI  (Z(J)-ZTC)/WID  l-OARGBI /(  DARGT- OARGB  ) B 

BF( J)-3AN8*0EXP(-(DEP-ZC Jl I/Q3EPI  B 

IF  (Ud.JI.LE.O.)  0(  l, J ) *BF( J)  B 

IF  (QCl, JI.EQ.O.)  OCltJI-BF(j)  B 

IF  (U( It JI .LE.O. » T ( 1 » J I -8  ( J I B 

03  120  1-2,  II  B 

U( I . J l-OBNDC  J)  B 

T(I,J)-8(J)  B 

Q(l,J)«BF(J|  B 

U(I,J>-0  B 

120  CONTINUE  B 

BVNAX-3  B 

00  130  J-2.JJ1  8 

BVNAX-ONAXlt  BVMAX,  ( C BC  J ♦ 1 » -B(  J » » /3Z  21  J 1)1  B 

A-DNAXICICUC  lll,J+ll-J(  Ill,J-ll)/DZl(J)/2)**2-(T(II,J+ll-T(II,J-  B 
11J>/021I J)/2,ZER0)  B 

130  Q(II,J!*BF(J!»A*RLAM*RlAN/COF  B 

BVNAX-3SQRTC BVMAX)  B 

Q(II,JJ>-Q(II,JJ1)  B 

0(11,11-0(11.2)  B 

00  1 AO  J-ltJJ  B 

00  140  1-2,111  B 

140  0(I,J)-0(11,JI  B 

DO  150  J-l.JJ  B 

D3  ISO  1-1,11  B 

ZU(  I,  JI-UI  I,  J)  B 

ZW(  I,  J)-W(  I,  J)  B 

ZTC I, JI-TI  I, J)  B 

Z0(I,JI>0( I,J>  B 

ADV( I , J 1-0.  B 

F 1 1 , J ) -0.  B 

Ed.JI-O.  B 

SU.JIO.  B 

A V2I I » J)  -0,  B 

P ( I , J 1-3  B 

PEI  I,  JI-PC  I,  J)  B 

150  continje  b 

C B 

C B 

C TINE  STEP  AN3  COUNTS  B 

C B 

OT-CFL/I  UN  AX  / DX1  (2  ) ♦BVMAX)  B 

TINE-TINE*G4P/UNAX  B 

NKN-T IN=/N0UTPT/DT/2*.5  B 

NKN-2-NXM  B 

NN-NXN*  '(O'JTPT  B 

DT-TIHE/NN  B 

NLF-2.0  B 

Ni.F-3  B 

NC-0  b 

nl-nlf  b 

NK-0  B 

WRITE  (6,160)  TIt€ , DT , NN  ,N  KN  ,NLF  B 

160  FORNAT  (//•  TOT  I HE  , 3T*,2FI2.4/*  NN ,NKH, NLF • , 3 16/ / I B 

TINE-3.3  B 

C B 

C B 

C INITIALISE  TINE  STEPPING  ARRAYS  B 

c e 

TP3-2.3  6 

03  1 T J I«1,M1  B 

OXP  ( I ) « 3T  / C OX  l C I *1 1 *DX  2 ( II  I 6 

OXN(  I)OT/(3Xl(I)*DX2l!  I)  B 

DXT(lt-3T/(2*0Xl(IH  B 

OXG(I)*2*OT/DX2(I)  B 

OXU( I l-3T/(4-0X2(l I)  B 

OXWI I l«3T/(4-0Xl(I | I B 

X*OR( I l-l-P0RC*0T*BVM4X*0EXP(-3S INH( DSINHIOMINl (6-6-RI I ( /3AP.TP3  8 

nm  b 

170  CONTINJE  B 

C B 


182 

183 

184 
1B5 
186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 
22S 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 
24  3 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 
25B 

259 

260 
261 


l 


25 


00  160  J-l.JJl 
02P(J)>)T/(DZ1(J+1I«0Z2(JI  I 
0ZM(JIOT/(DZl(J)*DZ2(  Jll 
0ZT(J>«7T/(2*DZ1(J>> 

DZG( J I - 2*D7/DZ2( J) 

ozai  ji»  jt/u-dzmji  ) 

OZWI  J)OT/(4«DZ2(JI  I 
140  CONTINUE 

C 

c 

C FOR  POISSON 

AAR-P I «? I/GAP/GAP/ 8 

A AZ >011(1(4/ 0X1(21 P-2-AAR, 4/OZlI 21  **2* AAR, 4/DZ 1 ( J J l ) --2+AAR > 

00  HO  1-2,111 
All l 1 -1.000/ (0X1 ( l I-DX2I  II  I 
A2( 1 1-1.000/(0X1  (I I-DX2 ( I- 1 » 1 
A3(  1 )« 1 1 ( I )*A2  ( I l- AAR 
HO  CONTINUE 

00  200  J-2.JJ1 

B1  ( J ) « 1 .000/ ( 0 Z1 ( J ) *DZ  2 ( J)  I 

B2( J) -1.000/ (DZ1 (J  >*0Z2( J-II ) 

B3(  J)»U(  JI+82IJ  HAAR 
200  CONTINUE 

C 
C 

CALL  R3»T 
C 
C 

C • *•«*«•  MARCHING  PTOCESS  BEGINS  **••*•* 

C 

c 

210  CONTINUE 

IF  { (L- NLF ) 260.220,260 
C 

c initialise  z f and  re-initialise  to  remove  parasites. 

C p IS  BACK  EXTRAPOLATED. 

c 

220  CONTINUE 

NL-0 

00  240  1-1,11 
00  230  J-2.JJ1 
PST-PEI  I,JI 
PE(J,J)»2*PST-P( I, J I 
P ( I , J ) »*ST 
230  CONTINUE 

00  240  J-l.JJ 
ZU( I, JI»U( I,J) 

ZWII, JI»WII,JI 
ZT(  I,  Jl-Tl  I,  J) 

Z0(I  , J)-0(  I,  J) 

240  CONTINJE 

CALL  HATCH  ( U, ZU ,W , ZW, T , ZT , 0, ZJI 
'JO  25  0 J-l.JJ 
00  250  I -1  , 1 I 

ZUt  I.  JI*(UI  I ,J  l*ZU(  I.  J)  1/2.0 
Z4(  I,  J*  »(H(  I , J l*Z4(  I,  J I 1/2.0 
ZT  ( I , Jl«( T ( I , J » ♦ ZT ( I ,J)  1/2. 

ZQ(  I,J)-(9( I , J ) + Z3 ( I.J) !/2. 

250  CONTINUE 

C 

260  CONTINJE 

c 

C 

CALL  MARCH  ( ZU.U.ZW.W.ZT.T.ZQ.Q) 

CALL  MARCH  (U.ZU.W.ZW.T.ZT .Q.ZOI 
TIME- T I NE»DT*2 . 

NL-NL 

NC-NC»2 

NK-NKAZ 

C 

C FORM  OUTPUT  *••••** 

c 

CALL  OUTPUT  (U.ZU.M.ZW.T.Q.F.E.S.ADV, R.X.Z.Y.DX l.DZl.RLAN.OT.TIN 
1E,GAP,0E»,I1 ,JJ,NC,NK,NKM.I2 , Jl, B, 0X2 , DZ2 , TT , T B, ZO , DSC ALE, NN , T I N.W 
2 IDFT.AJV, AV2.NC0PTSI 

c 

c 

C MODIFY  PLANT  OUTFLOW  TEMPERATURES. 

C 

T FAHR  ( 2J-T  INI  1 MTFAHRI  II 


B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

6 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 


262 

263 

264 

265 

266 
267 
266 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 
2B0 
281 
282 
283 
264 

285 

286 
287 
2B8 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 
326 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 


26 


• 

TFAHR (3)-TIN(2 )*TFAHR(4) 

6 

342 

00  280  J-l.JJ 

8 

343 

( 

IF  (LM  1, J) .LE.O. ) CO  TO  280 

8 

344 

A-l.0-10 

B 

345 

0-0 

B 

346 

00  270  K-1,4 

6 

34  7 

AA-FUNC  ((Z(J)”ZCENTIK))/'APL(K)  ) 

B 

348 

AU«AMPJ(K)/RPL(K) 

B 

349 

AV-DMAX l ( AU, ZERO) 

8 

350 

A-A-AV-AA 

8 

351 

r 

0*D*AV«( TFAHRt K )-T B > / ( T T-T 8 ) *44 

6 

352 

270 

CONTINUE 

B 

353 

T(1,J)0/A 

8 

354 

ZT( 1, J)-0/A 

8 

355 

280 

CONTINUE 

8 

356 

C 

B 

357 

c 

8 

358 

• 

IF  (NC-NN)  210,290,290 

B 

359 

i c 

c 

8 

360 

290 

CONTINUE 

B 

361 

c 

8 

362 

STOP 

6 

363 

* 

END 

FlWCT ION  FUNC  (A) 

C 

1 

r 

c 

»** 

FUNCTION  FOR  PLANT  INFLOWS  AND  OUTFLOWS 

C 

2 

I9PLICIT  REAL-81 A-H.0-2 ) 

C 

3 

8-1.00-APA 

c 

4 

IF  (8.LT.0.00)  8-0.00 

C 

5 

FUNC»B«8 

C 

6 

RETURN 

c 

7 

END 

« 

SUBROUTINE  MARCH  ( U , ZJ, W , 2 M, T, 2T , Q , 2 0 ) 

0 

1 

c • •• 

MARCH  ONE  HALF  TIME  STEP  — 

0 

2 

IMPLICIT  REAL-81 A-H.0-2) 

0 

3 

c 

0 

4 

COMMON  /N  l / OT.T  IME 

0 

5 

COMMON  /N2/  COF.RLAM.PI ,CAP,0:P,CSTRAT,CW 

0 

6 

COMMON  /NO/  II,JJ,II1,JJ1, II2.JJ2 

0 

7 

€ 

COMMON  /NT / A1 ,A2,A3,B1 ,82 ,83,WI ,WK 

0 

8 

COMMON  /N1 2/  F , E ,A DV«4  V2, 5 , P,Pc 

D 

9 

COMMON  /N1 3/  ALP,BET,CAM,t)EN,H,S0,2S,S3l,S02,2Sl,2S2,ZFTU, 2FTM.2 

D 

10 

1FQU.ZFQM,  <FTU»  XF  TM ,XF QU.XF  ON ,XPOR 

0 

11 

COMMON  /N14/  OX1,OX2,OXP,OX9,OXT,OXC,OX3,OXU,OXW,321,D22,32P,DZM 

0 

12 

1. 

02T,02C,02D,02U,DZW,2BL,2BR,2WL,2WR 

0 

13 

COMMON  /N15/  R , X,2 • Y 

0 

14 

COMMON  /N16/  B , BF,  UBNO,  F AC 

C 

15 

C 

z 

0 

16 

REAL  abs.sort 

0 

17 

real  T,2T,Q,ZQ,F,E,ADV,AV2 

0 

18 

REAL  U,2U,W,2W,S.P,PE 

0 

19 

c 

0 

20 

c 

0 

21 

DIMENSION  U( 12  ,401 , 2UI12.40I,  W(12,40l,  2W(12,40),  S(12,40),  P( 

n 

22 

112,90),  PE ( 12, 40 ) 

0 

23 

« 

DIMENSION  T ( 12 ,40)  • 2TI  12.40),  0(12,40),  20(  12,40),  FU2.40),  E( 

0 

24 

112,40),  4V2( 12,40) , ADV(12.4l)) 

0 

25 

DIMENSION  R(  25  ) , X(25),  2(  90),  V(40) 

0 

26 

DIMENSION  0X1(25),  0X2(25),  0XPI2SI,  OXM(25),  DXTI25),  DXS(25), 

0 

27 

10X0125),  0XU(25>,  DXW( 25) 

D 

28 

DIMENSION  021(40),  022(47).  02P(40),  02M(40),  D2T(40),  023(40), 

0 

29 

1020(40),  020(40) , 02W(40),  Z3U40),  Z8R(40),  2WL(40),  2HR(40) 

0 

30 

DIMENSION  Al(25>,  42(25  ),  43(  25),  81  (40),  B2I40),  B3I40).  MK25) 

0 

31 

• 

1 , 

WKI40) 

0 

32 

• 

DIMENSION  8(40),  FACI43),  BF(43I 

D 

33 

DIMENSION  UBNO (40) 

0 

34 

OIMENStQN  AL  P ( 40 ) , BETI40),  349(40),  0EN(40),  H(4DI,  S0I43),  2SI 

0 

35 

1401,  SO 1( 401 , SO  2( 40 ) , 2 S 1 ( 4 3) • 2S2(40),  XFTU(40),  XFTMC40I,  XFQU( 

0 

36 

240),  XF  39  ( 40 ) , 2FTUI40),  2FT9(43),  2F0U(40),  2F09(40),  XP0R(40> 

0 

37 

c 

0 

38 

c 

0 

39 

• 

DT2-2-0T 

0 

40 

JJ1-JJ-1 

0 

41 

c 

0 

42 

CALL  MR  ( Q, * UBEF ' ) 

0 

43 

CALL  MR  (T.'TBEF*) 

0 

44 

1 

OJ  20  t-l,II 

0 

45 

t 

— --  z^r-z- 

27 

J 

f 


T 


C 


e 


c 


c 


c 


€ 


c 


c 


03  10  J-2.JJ1 

0 

46 

PST -PEI I,J) 

0 

47 

PEI  It  J)«PI  l« 41 

D 

48 

Pll»J)»2*PII»J  l-PST 

0 

49 

FI  I,  J)-RLAM*SQRTI3U.J)  )/l  1.«-RLAN*RLAM*ITI  I,  J«-l  l-TI  I , J-ll ) /0Z1I J 

0 

50 

1)/2.*CSTRAT/QI t,J)  1 

0 

51 

c 

02/)  »l  J,IIWISBA4|)J,I(UtSBA< 

D 

52 

c 

D 

53 

IF  (FII,J).LT.O.)  F(I,J)-0. 

0 

54 

10 

CONTINUE 

0 

55 

FII,l)>FII,2) 

0 

56 

FII,JJ)-FII,JJ1) 

0 

57 

23 

CONTINUE 

0 

58 

CALL  WR  I F t • F •) 

0 

59 

C 

0 

60 

C ZBL 

ANO  ZB*  ARE  1 FOR  INFLOW,  0 FOR  SLOW  OUTFLOW , -l  FOR  FAST  OUTFLOW. 

0 

61 

C ZWL 

AND  ZWR  ARE  1 FOR  INFLOW,  0 FOR  SLOW  OUTFLOW,  -1  FOR  FAST  OJTFLOW. 

0 

62 

AA-  3 

0 

63 

Cl-0 

0 

64 

0-0 

0 

65 

ZEAO-O. 

0 

66 

ONE-l. 

0 

67 

00  30  .1-1, JJ1 

0 

68 

A-Ull.J) 

0 

69 

C-FI2,  J 1/0X1  12  )♦  A/ 2. 5/ 100 

D 

70 

ZBLIJI-0.00 

0 

71 

IF  (A. ST. ZERO)  ZBL I J )-3NE 

0 

72 

IF  1C. LT. ZERO)  ZBLIJI— ONE 

0 

73 

A-A-AA 

0 

74 

A A- A- A A 

0 

75 

C-C-Cl 

0 

76 

Cl-C-Cl 

0 

77 

ZWLI J) >3.00 

0 

78 

IF  IA. ST. ZERO)  Z WL 1 J ) -ONE 

0 

79 

IF  1C  .LT  .ZERO)  ZWLIJ)  — ONE 

0 

80 

A-UII 11, J) 

0 

81 

C-FI  I Il,J)/3XltIU)-A/2.5/  100 

0 

82 

Z 3R I J ) - J.DO 

0 

83 

IF  14. LT. ZERO)  ZBR 1 J)-3NE 

D 

84 

IF  1C. LT. ZERO)  ZBR 1 Jl --ONE 

0 

85 

A- AAA A 

0 

86 

aa-a-aa 

0 

87 

C»C*C1 

0 

88 

Cl-C-Cl 

0 

89 

ZWRI Ji-0.00 

0 

90 

IF  IA.LT. ZERO)  Z WR 1 J ) -ONE 

0 

91 

IF  (C.LT.ZER3)  ZWRIJ)  — ONE 

0 

92 

30 

CONTINUE 

0 

93 

C 

D 

94 

C 

0 

95 

IF  ITME.E3.0.)  CALL  BOUND  I ZJ  , ZW,  ZO,  ZT , Z BL  , ZBR  , ZWL  • ZWR  , B,  BF  , 1 1 , 

0 

96 

lJJ.OZl) 

0 

97 

IF  IT  I NE.EQ.O.  ) CALL  BOUNO  I U,  W, 0, T , ZBL  , Z3R  , ZWL , ZWR,  B,  BF , I I , JJ ,0 

D 

98 

1Z1) 

0 

99 

c 

0 

100 

c 

D 

101 

CALL  FCSTD  1 U,W,T,ZT,3,Zg,F,A0V, AV 2,  ALP , BET , CAM , DEN,  H,  SO l,  S02,  ZS 

0 

102 

ll,ZS2,ZFTU,ZFTM,ZFQU.ZFQH,XF TU.XFTh.XF OU, XF3M, 0X0, OXP, DXM,DXT,E,S, 

0 

103 

2ZBL,ZB*,3ZD,OZP,OZM,OZT,3XI,OX2,OZ1,OZ2,FAC,BF.B, I I, JJ I 

0 

104 

CALL  FCSTU  1 U, ZU ,W,F ,E , AOV , ALP ,BET .CAM, OEN ,H, SO, ZS,0XP,0 ZP , ZBL ,Z 

0 

105 

1BR,0XN,DZW,0XU,0ZU,I I, JJ.P.0X2) 

0 

106 

CALL  FCSTW  1 W,ZW,U,T,F,E,AOV,ALP,BET,CAM,OEN,H,SO,  ZS ,OXP,OZP,Q,X 

D 

107 

1P0R,ZWL,ZWR,DXN,0ZM, 0XW.3ZW,  1 1 , J J , P, 0Z2) 

0 

108 

CALL  WR  1 ZT, • T • 1 

0 

109 

CALL  WR  IZQ, ’Q  • ) 

0 

110 

CALL  WR  IZU.'Ul  * ) 

0 

111 

CALL  WR  IZW.’Wl  •) 

0 

112 

c 

0 

113 

SIGN— l. 

0 

114 

CALL  TSTEP  1 ZU, ZW, P, OXC.OZC, II, JJ, SI CNI 

0 

115 

c 

D 

116 

00  40  1-2,111 

0 

117 

30  43  J-2, JJ1 

0 

118 

SII.J)  — 1 ZU  1 1 , J ) — ZD  I I — 1 • J ) ) /DX1  (Il'PIZWI  I.JI—ZWII,  J—  1)  I/0Z1I  J) 

D 

119 

S(I,J)-SII,J)/DT2 

0 

120 

40 

CONTINUE 

0 

121 

c 

D 

122 

CALL  PUSS  IP,S,E,A1,A2,A3,B1,B2,B3,WI,WK,II, JJ) 

0 

123 

c 

0 

124 

c 

0 

125 

€• 


28 


o o o o o 


c 


f 


SIGN-1. 

D 

126 

CALL  T STEP  I 2U  ,2M,  P ,OXG,OZ  5,  II, JJ, SIGN) 

0 

127 

CALL  BOUND  ( 2U  ,ZW,  23 ,2T  , ZBL,  2BR,  ZML.  ZMR,  B.  BF,  1 1 , JJ , 02 1 1 

0 

128 

00  50  1-2,111 

0 

129 

00  50  J-2.JJ1 

D 

130 

S(I«J)>(ZU(I,J)-ZU(I-1,JI) /0X1 (l)-(ZM(I,J)-ZM(t,J-l))/DZl(J) 

0 

131 

S(I,J)-S( I.JI/0T2 

0 

132 

50 

C0NTINJ5 

D 

133 

C 

0 

13A 

c 

0 

135 

CALL  MR  (P.'P  •) 

D 

136 

CALL  MR  (ZU, *UF  • 1 

0 

137 

CALL  MR  (ZH.'MF  •) 

0 

13B 

RETURN 

0 

139 

END 


c 


r 


c 


c 


c 


c 


« 


SUBROUTINE  FCSTD  ( Ui  M,  T , ZT  , <J,ZQ,  F , ADV,  AV2 , ALP . BET,  GAM.OEN,  H,  SOI , E 
1S02«ZS1«ZS2,ZFTJ,ZFTM,ZF3U»ZF<1M,  XFTU»XFTM,XF3U»XFQM, DX0»DXP,9XM,DX  E 


2T,E. S , Z8L. 29R.  DID,  02 P.OZM,  02  T,  0X1,0X2,021,02 2 ,FAC  , BF  , B,  ML.  NL  I E 

C •••  FORECAST  T ANO  0 • ••  E 

IMPLICIT  REAL-81 A-H.0-2I  E 

C E 

COMMON  /HU  DT.T1NE  5 

COMMON  /N2/  COF,RLAM,PI ,GAP,DEP ,CSTR AT ,CM  E 

C E 

REAL  ATS, SORT  E 

REAL  T,2T,Q,2Q«F,ADV,AV2,J,j,E,S  E 

C E 

DIMENSION  U(NL ,NL> , W ( ML ,NL  ) , =(ML,NLI,  S(ML.NL)  E 


DIMENSION  T(NL.NL),  3<ML,NL>,  F(NL.NL),  ALPCNL),  BET  I NL ) , CAMJNL  E 
1),  DEN(NL),  HIND,  SOl(NL),  S02INLI,  ZSHNL),  ZS2INL),  ZFTUI  NL  1 , 2 E 
2FTMINLI,  2FQUINL),  2FQMI NL I , XFTJ(NL),  XFTMI NL I , XF3UINLI.  XF3MINL  E 
31,  OXOIML),  OXPIMLI,  OXMINLI,  OXT(ML),  OZDINL).  02PINL),  02MINL),  E 
ADZTINL),  0X1 (ML ) , 0X2IMLI,  0 2 1 1 NL ) , 022INLI,  FACINL),  BF I NL ) , AOVI  E 
SML.NU,  XV  2 1 ML  ,NL ) * 2T  ( ML  ,NL  ) , ZOIML.NU,  23LINL),  ZB»(NLI,  BIND  E 

E 

E 


2ER0-0  E 

TMO-2  E 

II- ML  E 

JJ-NL  E 

III- ll-l  E 

II2-II-2  E 

JJl-JJ-l  E 

JJ2-JJ-2  E 

E 
E 

2 IMPLICIT  E 

DO  10  J*2 , JJ 1 E 

XFTUCJ)-(F(2,J)-F(  l,J))/2.J0J-(ZT(2,J)*ZT(  1 , J ) ) E 

XFOUI JI-(F(2,J|*FC  l.JII/2.  300*  C 201  2,  J I -231  1,J1)  E 

10  CONTINUE  E 

C E 

00  70  1*2,111  E 

C E 

DZO(l)»(FI  I,  2) -Fit,  1)1/2.000  E 

00  20  J*2, JJ1  E 

0 ZQ  < J I * ( F < I , J ♦ 1 1 *F  ( I , J ) 1/2,003  E 

23  CONTINUE  E 

00  33  J*2, JJ1  E 

FA*(F( t»l,J> »F ( I, J ) 1/2.003  E 

XFTM(JI*XFTU(J  I E 

XF3MI Jl «XFOU( J ) E 

XFTU(JI-FA-(ZTIl-l,JI-ZT(I ,J)>  E 

XFQUI  JI-FA-<ZQ(I-l,JI-ZO((,J>)  E 

30  CONTINUE  E 

C E 

SAMI 1 1 *-020(1 1 *OZM III  E 

00  AO  J*2,JJl  E 

AA»CJ(  t,JI-J(  I-1,J  I l/OXK  I » E 

AB*(  J( I,  J»1)*UU-1  ,J*1I-J( l, J-ll-JI I-l.J-l II/02KJ l-IMM-l ,J|*W<  E 
1 I-l.J-l) -Ml  I-1,J)-W<  I-l.J-l)  )/OXl(l>  E 

AOVI I , J)*OXT(I )* (U 1 1 , J I *T I l-l,J>-JU-l.JI*T(I-l,J>l-DZT(J)*<W(I,  E 
1JI-TU,  J-ll-WI  l,J-l)*TII,  J-D)  E 

C E 

C UPSTREAM  FINITE  DIFFERENCE  REPRESENTATION  FOR  3.  E 

UO*UII  .Jl-Ut  1-1, Jl  E 

MQ-MI  I , Jl -Ml  I , J—  1 1 E 

10-1  E 


1 

2 

3 

A 

5 

6 
7 
B 
9 

10 

11 

12 

13 

1A 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 

25 

26 
27 
26 

29 

30 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 
Al 
A2 
A3 
AA 
A5 
A6 
A7 
A6 
A9 

50 

51 

52 

53 
5 A 

55 

56 

57 
SB 

59 

60 
61 
62 
63 


( 


29 


r 


c 


c 


c 


t 


« 


t 


JO*  l E 

IF  (U3.LT.0)  10—1  E 

IF  (PO.LT.O)  JO— l E 

: AV2II  • JI-OXU II*(U(  I,J)*0( t*l, JI-JI 1-1, J)*Q(I-1. J) )*D2TI JI*(WI I,  f 

C IJ)*Q(1,  JUI-PI  I,  J-l)*Q(t,  Mil  E 

PB*(P((,JI«W(t,J-l))«.SD3  E 

4V2( I, UI-UQ*OXT( 1 ) *1 0*  < 0 < I .JI-OC I-IO.J) >*PQ*0ZT< Jl  • JQ* ( 0 ( I , J l-OI  E 

11,  J-JOI )-3T*F(I, Jl*( (a.*4A*A4»A5*A8/16.-(TU,J+L)-T(  I , J- 1 1 I / 2. /021  E 

2 I Jl l-QI I,J>*C0F/ALA.N/ALAN>-P8*PBPCP*F( I,JI*3T*(T( I ,J»1)-T( I, J-l) 1/  E 
32./0Z1I JI/QI I.J)  E 

*4*0.03  E 

*8*0.00  E 

*»OT*F(  I,  J )*({*.**  A*AA»AB*  *8/ 16.  ♦(  Ul,  J*l)-T(  I , J-l )) /2./D2  II  Jl  I / E 

1 Q(  I,JI*C3F/RLAM/RLAM)*U0*I3*3XT(I)*MQ*JQ*DZTIJ)  E 

ALPI  Jl— 02D(J-ll*0ZP(J-l)  E 

04*41  J I — DZD(  J)  *D2NI J ) E 

8ETI J t *1 .000-0 AMI J l-ALP  I J)  E 

E( I , J ) *4  E 

SOU  J)*ZTt  I,  JI+OXMI I |*XFTJ  C JI-OXPI 1-1 ) •XFTNI J I - AOV C I , J I E 

S02 I J I *ZQ( I,  J) *DXH( I )*XF JU( J)-DXP( I- II *XFQN( J I - AV2 ( I » J)*2* A*Q( I,  E 
1J>  E 

AO  CONTINUE  E 

c E 

BETI2I*BET(2)*ALPI2)  E 

BETIJJII-3ETI JJ1  }»GAMI JJ1)  E 

C E 

CALL  T X 1 01  C ALP, BET, GAN, SOI, 2S1.JJ1  E 

03  50  J*2, JJ1  E 

BET(J)*8ET(J)*E(  !,  J I *TPO  E 

50  CONTINUE  E 

CALL  MIDI  I ALP,BET,GAN,S3  2,2S2,JJI  E 

C E 

00  60  J*2, JJ1  E 

ZT( I, JI-2S1IJI  E 

S 1 1 , J I *2Q(  I, J I E 

231 1 , JI*0NAX  II  ZS2I Jl ,2J( I, JI/TP3)  E 

63  CONTINUE  E 

C E 

2TC 1 , 1 1 *2  T I I ,2)  E 

2TI  I.J  JI-2TI  I,  JJ1I  E 

23(1.11*2311,21  E 

231  I.JJI-20I  I.JJll  E 

73  CONTINJE  E 

C E 

C E 

CALL  PA  C 2T,  • T INT*  I E 

CALL  PA  (20, *3  INT* ) E 

CALL  6)  (E.ML.NL)  E 

CALL  PA  (E  ,•  E-A  'I  E 

C E 

C X IMPLICIT  E 

00  83  1*2,111  E 

2FTUI  I l-CFU  ,2  )*F(  1,1)  1/2.  303*(2T(I,2I*2T(  1,1)1  E 

2F3UC I)-(F(I ,2)*F(  1, 1)1/ 2. 300* (20(  1, 21-231  1,1)  1 E 

83  CONTINUE  E 

C E 

C E 

00  1 AO  J*2 , J Jl  E 

c c 

C E 

00  93  1*2,111  E 

2FTNI I l*2FTJ( I ) E 

ZFJNI I )*2F0U( I I E 

FA*IF(I,J*l)*FII,JI)/2.000  E 

2FTU(II*FA*(2TII,J*1I-2T(I,J1I  E 

2F3'JI  I)-FA*(20(I,J«1)-Z0(I,JII  E 

93  CONTINUE  E 

C E 

03  100  1*1,111  E 

0X0 ( I)*(F( l*l,J)*FI I.J) 1/2.033  E 

130  CONTINUE  E 

C E 

03  110  1*2,111  E 

GAN 1 1 ) — 0X0( 1 1 *OXN ( I ) E 

ALP1II—  0X01 I-ll*OXPII-ll  E 

3cT(II*l. 300-G AN (I)-ALPII)  E 

SOU! )*2T(I,JI*02N(JI*2FTU(I)-02PIJ-l)*2FTN(l)-A0VII,JI  E 

502(1 l*20( I, JI*02N(J)*’  >U(I)-02P( J-I)*2F0N( I)-AV2(I ,J)*E( I, J)*(  E 

12. *3( I • JI-SI  I, J)  ) E 

110  C3NTINUE  E 


*6 

65 

66 

67 

68 

69 

70 

71 

72 

73 
76 

75 

76 

77 

78 

79 

80 
81 
82 
83 
86 

85 

86 

87 

88 

89 

90 

91 

92 

93 
96 

95 

96 

97 

98 

99 
100 
101 
102 
103 
106 

105 

106 

107 

108 

109 

110 
111 
112 
113 
116 

115 

116 

117 

118 

119 

120 
121 
122 
123 
126 

125 

126 

127 

128 

129 

130 

131 

132 

133 
136 

135 

136 

137 

138 

139 
160 
161 
162 
163 


I 


c 


30 


w 


r 


r 


c 


€ 


c 

C LEFT 


C 

C RIGHT 


C 


120 

C 


130 

1*0 


IZBLIJI.GT.ZER3) 
(Z  dL  < J ) . GT  . ZER3  I 
(ZlLlJ).Ea. ZERO) 
CZSLIJI.LT. ZERO) 
IZSLCJC.LT. ZERO! 


S3l(2)*S0l(2)-ALP(2l*ZT(l,J) 
S32(2)-S32(2)-ALP(2>*ZQt l.J) 
BET (2) "SET ( 2 ) »ALPI 2 ) 

BET  I 21 -SET(2)*ALP(2)«2.00 
GANI2)-GAM|2|-AlPC2I 


IF  IZ5RIJI.GT.ZER3) 
IF  CZBRC JI.GT.ZEROI 
IF  IZ8RCJI.EQ.ZER3  > 
IF  CZSRC Jt.LT.ZER3) 
IF  IZSRCJI.LT. ZERO) 


S01CII l)*S3l( Ill)-GAMC  1 1 1 ) *B C J I 
S32CII  1I-S32I  Ill)-GAN<  1 1 1 1 *BF(  j) 
BET  C II  l)-SET(Ill)  ♦GAMC  1 1 1 1 
SET  CII 1I-3ETCII1) ♦GAM) I ID *2.30 
ALP  C II  D-ALPCIID-GANI  DU 


CALL  TRIQl  (ALP, BET, GAM, S01.ZS1.II) 
33  120  1-2. Ill 
BETIll-JETCIDEI  I.J) 

CONTINUE 

CALL  TRI01  C ALP. BET. GAN, S32.ZS2.lt  I 


03  130  1*2.111 
ZTCI.JI-ZSIC I) 

ZERO-OCI ,J)/2 
EPS-l.3-20 

Z3C  I,  J)«0MAXIC ZS2C  D.ZER3,  EPS) 

CONTINUE 

CONT I N IE 

RETURN 

E NO 


E 144 
E 143 
E 14* 
E 14T 
E 148 
E 145 
E 150 
E 151 
E 152 
E 153 
E 154 
E 155 
E 156 
E 157 
E 158 
E 159 
E 160 
E 161 
E 162 
E 163 
E 164 
E 165 
E 166 
E 16T 
E lcB 
E 169 
F 170 
E 171 
E 172 


- 


r 


c 


c ' 


c 


c 


c 


SJSRJUTINE  FCSTU  CU.ZU.M.F ,E , AON , ALP , BET , G AN, DEN ,H , S3, ZS , 3XP , 0 ZP 

1 . ZBL.ZBR, JXM, OZN.OXU.OZU.NL, NL.P, 0X2) 

C •••  FORECAST  U ••• 

INPLICCT  REAL-BCA-H.O-Z ) 

C 

COMMON  /Nl/  DT .TIME 
C 

c 

REAL  F.E.AOV 
REAL  'J.ZU.M.P 
C 

0IMENSI3NUCML.nl),  ZUCML.NL).  <C(ML,NL),  FIML.NLI,  EC-L.NL),  aov 
1 CML.NL ) » ALPCNL).  BETCNL),  GARCNLI,  OEN(NL),  HCNL),  SOCNLI.  ZSCNL) 

2,  OXP(ML),  3ZPCNL).  OXMIML),  OZMCNL),  OXJ(ML),  DZU (NL ) , ZBLCNLI,  Z 
3BRCNL).  * ( ML  .NL  I . 0X2  ( ML  I 

OIMECSIJN  ZLC200I 

c 

c 

II-ML 
I II -I  1-1 
112*1 1-2 
JJ-NL 
JJ1-JJ-1 
J J2-JJ-2 
C 

00  10  J«1,JJ1 
30  10  1*1.111 

El t.J)*IF( l»l.J*l)*F(I,J.l).F( I»1,J)»F( I , J ) 1/4.0 
10  CONTINUE 

C 

c 

c Z IMPLICIT 

00  20  J-2.JJ1 
ZLI JI-ZJI1.JI 
20  CONTINUE 

C 

DO  SO  1*2.111 
00  33  J-2.JJ1 
GAMI J I — OZNI JI-EII.J) 

ALPCJ)*- OZPCJ-ll-EC I.J-ll 
BET  I J ) *1  .000-GAM  I J I- ALP  I J) 

OKP*OXP| I )*F 1 1 41 ,J I 
DKNOX4I  I)*FCI  ,J) 

TJX*  IUI  Dl.JDUC  I,  JI)*UC  IM.JI-CUC  I . J I » Jl  I - 1,  J ) l-UC  t-l.JI 
TUZ*l<CII»l.JI*MC  I,  J ) ) *U  ( I , J*  l ) -(  M(  I ♦ 1 » J- 1 1 ♦*  ( I « J-l  ) ) *U ( 1 , J -l  I 
AOVI I , JI-OXUI I »-TUX*OZU(  JI*TUZOT*(P(  1H.JI-PCI,  J)  )/DX2t D-16.-0 
1XUI I )«0Z  JCJI/OT- IE  II .J )*IJI I »1,J)-M( I, J ) I - EC t,J-l)*(W( !♦!, J-ll-MC I 
2.J-1II) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


c 


31 


r 


€ 


€ 


r 


c 


< 


€ 


FX«DRP-ZU(I-1,JI-0KM*ZL(J1“(D<P-0KN)*ZJ(  It  J» 

S3( J)«ZJII ,J)-AOVI t,J)+FX*2. 

ZL(J)-2)(I,J) 

33  CONTINJE 

C 

BET(2)-3ET(2)»ALPI2) 

dET(JJL)-BET(JJl)-GAM(JJl) 

C 

CALL  TU01  (ALPt8ET,CAH.S3t2S.JJI 
C 

00  40  J-2, JJ 1 
20(1, J)-ZS(J) 

40  CONTI  HUE 

2JI  I,  U-Z'J(  1,2) 

ZUC I.J JI-ZUI I.JJI) 

53  CONTINUE 

C 

C 

CALL  AT  (ZUt  'UINT*  ) 

C X IMPLICIT 

00  60  (-2,111 
ZLI 1 1 -231 1,1) 

60  CONTINUE 

C 

00  93  J-2.JJ1 
33  73  1-2,111 
ALPI  1 )—  OXMt ll-FCI  . J)*2. 

5AM(  l|.-DXP(  I l*F(I-l,J>  -2. 

BET  I I )-l  .ODO-ALPd  l-CAHC  I) 

OKP-OZHI J) *E ( I , J> 

OKN-32  °( J- 1 )*E(I,J-l) 

F2*3KP*2U(  I,  J-ll-DKM-ZL  ( ll-(0<P-0<M)*ZU(  I,  J) 

SOI I)-ZJ(I ,J)-AOV( I.JI-FZ 
2LID-2JII.J) 

70  CONTINUE 

S7I  2) -SOI  2 )>ALP(  2)  *ZU(  l,  J) 

IF  IZ3UJI.CE.0.)  ALPII  Il)-ALP(Ill)  ♦GAH(  1 1 1 1 
IF  (ZBd J).LT.O.)  ALP(II1)-ALP(III ) -GAM (111) 

IF  (ZBR(J).LT.O.)  BETIIIl)-aET(IIl) ♦GAM) III) -2. 00 
C 

CALL  TRID1  ( ALP,BET,GAH,S3,ZS,II) 

C 

30  93  1-2,111 
201 ItJ)-ZS(I) 

93  CONTINJE 

>0  CONTINUE 

C 

c 

RETURN 

END 
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5 

6 
7 

a 

9 

0 

1 

2 
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SUBROUTINE  FCSTW  ( H, 14, J.T ,F.E . AON , ALP, BFT , GAM, DEN, H, SO, 2S ,DXP,0 
1ZP. J.XPOR.ZHLt ZHR.OKM.OZM.OX J, D2J,ML,NL,P,022) 

C FORECAST  4 

IMPLICIT  REAL-81 A-H.0-2 I 
C 

COMMON  /Nl/  DT .TIME 

COMMON  /N2/  C3F,RLAM,PI,GAP,0EP,CSTRAT,:h 

c 

REAL  U,P,M,ZM,Q 
REAL  F.E.T.AOV 
C 

DIMENSION  QIML.NL) 

01 ME NS  I ON  M(ML.NL),  24 1 ML, NL I,  J(ML.NL),  T(ML.NL),  FIML.NL),  E(M 
1 L , NL ) , A JV(  ML,  NL  ),  ALP  (NL  I » BET(NL),  GAM(NL),  OEN(NL),  HINLl  , SO(N 
2 L ) , 2S( NL) , OXP(ML),  OZP(NL),  OXM(ML),  OZM(NL),  OXM(ML),  QZrf(NL), 
3ZWLINL),  ZMRINL),  P(ML,NLI,  022INLI,  XPORIML) 

DIMENSION  ZL 12001 
C 

I I-  ML 

III- II-1 
1 12 - 1 1-2 
JJ-NL 
JJi-JJ-1 
JJ2-JJ-2 

C 

OTH-OT/l.ODO 

C 

c 
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cinpoooooopnonocicicincipnpoopopo 


c z implicit 

DO  ID  J-2.JJ2 
ZKJI-ZUll.J) 

10  CONTINUE 

C 

DO  40  1-2,  III 
DO  20  J-2, JJ2 

GANIJI—DZPI  JI*F(I , Jdl  *2.00 
ALPIJI— DZM(  JI*F(l,JI*2.00 
BET ( Jl-l.  ODO-ALP  (JI-GANI  Jl 
OKP-OXMI II*EII,JI 
OKN-OXPI 1-11*6(1-1 .J) 

FX-0KP*Z4I  I*1»JI*0KM*ZL(  J) -( OXP* DKNI *Z4(  I.  J) 

TUX«(U(l ,J*ll*U( I. Jll*4( 1*1. JJ-IUI l-l.J*ll*U( I-1.JII*W(  I-l.J) 

T UZ-I  Ml  I ,J*ll*M(  I.J)l*U<I.j4-l)-(Md.J)»M(I,J-lll*M(l,J-l) 

AOV< I , J 1-0X4  < I 1*TUX*DZW (JI*TUZ-OTH*(T( I»J»1I*TI l , J I - T ( 1 1 , J ♦ 1 1 -T I 
1 II. J) >»OT*IP(l ,J*1)-P( I. J I I/DZ2I J I -16. *0X41 I)*3Zm(J)/0T*(6(I . J I* (U 
2 1 1,  Jdl- J(  I.JII-EI  I-1,JI*(U(  I-l.J*- 1 »-U*  I-l.J  III 

SO(JI-ZM( I.JI-AOtfl l,JI*FX-C4*(FI I, J)»F( I, Jdl  >•(  T(  I, Jdl-Td.JII 
l/dl  I , Jl *9( I ,J*l )I*0T/DZ2( JI*ZW( I, J) 

ZL( JI-Z4II.JI 
20  CONTINUE 

C 

CALL  MIDI  ( ALP. BET, GAN, SO, ZS,  JJ1I 
C 

DO  10  J»2» JJ2 
Z4(  I,  JI-ZSIJ  >*XPOAI II 
30  CONTINUE 

40  CONTINUE 

C 
C 

CALL  MX  (ZU,  ‘MINT*  I 
C X IMPLICIT 

DO  SO  1-2, 111 
ZL( 1 1 ■ 3*0 
SO  CONTINUE 

C 

00  80  J-2.JJ2 
00  60  1-2,111 
GAN(  II—  OXMI  I)*E(t.J) 

ALPIII—  DXP(I-ll*E(  I-l.JI 

BET(i I -1.003-ALP( I I -GAM ( 1 1 *C4* (F ( I . J I *F ( I , Jd 1 I * ( T ( I , Jd )- T ( I , J l 
1 l/(Q( 1 , J l+Q( I > J* 1 ) l*OT/OZ  2 1 J I 
OKP-O  Z°  ( Jl  *F  ( I , J*l  I 
OKN-OZ  N ( J I *F  ( I ,J| 

FZ»OKP*ZH( I, J»1I*DKN*ZL ( Il-(0<P»0KMI *Zrf(  I , J ) 

SOI II-ZMII  .JI-ADVI  I » J I ♦ F Z* 2 . DO 
ZL(  ll-ZMI  I.JI 
60  CONTINUE 

c 

c LEFT 

IF  (ZMLIJI.GT. 0.031  S0(2l-S0m-ALP(2l*>m,JI 
IF  (ZMLI Jl .EQ. 0.001  BETI2I-BET(2I»ALP(2I 
IF  (Z4LIJI.lt. 0.001  BET 121-8 ET(2)+ALP(2I*2. DO 
IF  (ZML(  JI.LT.O.DOI  GAN  ( 21 -G  A 1(2  I - AL  PI  2 I 
C 

C RIGHT 

IF  UMRIJI.EQ.O.DOI  BET(IIll-BET(ltll*GAM(  (III 
IF  (ZMAUI.LT. 0.001  BETIII  1) -3  ETdIll+GAM(  1111*2.00 
IF  (ZUXIJI.LT. 0.001  ALPI1I1I-ALPIII1 1 -GAM ( 1 1 1 1 
C 

CALL  TX1D1  I ALP. BET, GAM.SO.ZS.il  I 
C 

DO  70  1*2, II 1 
ZUI 1, JI-ZSI I l*XPOR( 1 1 
70  CONTINUE 

SO  CONTINJE 

C 

RETURN 

END 


SUBROUTINE  TRI 01  ( A.B.C . S, Z.NI 
C •••  TRI 01  AGONAL  EQUATION  SOLVER  ••• 

C 

C SOLVES  TXIOI AGONAL  SYSTEM,  1-2  TO  N-l. 

C A.B.C  JNCHAN GEO.  AMS  S DESTROYED.  SJLJTION  IN  Z. 

c 

IMPLICIT  REAL*8( A-H.O-Z I 

DIMENSION  AINI,  B(NI,  CINI,  SINI,  UNI 


G 85 
G 86 
G 87 
G 88 
G 89 
G 90 
G 91 
G 92 
G 93 
G 94 
G 95 
G 96 
G 97 


H 

H 

M 

H 

M 

H 

H 

H 


33 


NH-N-l  X 

0-BI2I  M 

SI2MSI2I/0  H 

OH  10  1-3, NM  X 

ztt-n-ci  i-ii/o  x 

0- 3<ll-A<lMM  l-ll  X 

10  SIIMISIII-All  MSI  l-lll/D  X 

Z(NMMSINN)  N 

00  20  J“A«N  H 

1- NA2-J  X 

20  tm*siu*i(  imzimii  x 

RETURN  X 

ENO 


c 


€ 


C 


< 


C 


SUB* JUT  I ME  * OP  T 

C •••  INITIALIZE  POISSON  SOLVER  ••• 

MPLICIT  AEAL*S( A-H.P-ZI 
C 
C 

COMMON  /NM/  AAR ,AAZ, KAPPA 
COMMON  /NS/  RRI6Ai,LL 
0 ME  NS  I ON  ALPIB),  BETIBIt  OISJ) 

C 

c 

ALPI1MAAR 
5ETIIM4AZ 
RA-BETI ll/ALPI 11 
00  10  Ml, KAPPA 
ALPI  MIMIALPI  IMBETIIIMR.50J 
BET  I Ml  Ml  AL  PI  IMBETII  IIP.SOO 
10  CONTINUE 

al>alpikappa»d 
36-BET IKAPPA*! ) 

c 

Z NEKTON'S  ITERATION  FOR  R,  FOR  EXTRAPOLATED  INITIAL  APPROXIMATION 
C 

RK>OM 1 N l I 2*AL, (AL*6E )**.SOOI 

OJ  20  J>1.10 

AK*I 36-RKI  / ( 3E*RK) 

A-IBE-RKI/IBEpRK) 

BK-IRK- ALI/(RK«ALI 
3-1  l-AM*  1.500/ A*3-1/RK 
A>-2*JE/IB6»AK|»2 
BD-2**L/IAL»RKM*2 

F0-83/3K/3K-I 1-AK) ••.503*1 M4K/2 >/AK/AK*AD 
CM  R-l  AL»9SM».530  IF  NOT  EX T* AP 1L AT INO 
20  RK-RK-I/FO 

S-l/l  1/3K-3I 
AS-AK 

c 

C FINO  ALL  KK  A VALUES 

C RRIKI  ARE  THE  ITERATION  PARAMETERS 

C 

LL«2«*<4PPA 

RRtlMRK 

Kl-LL 

00  AO  MI, KAPPA 
J-KAPPA+l-I 
AL-ALPIJI 
SE-BETI Jl 
K2-K1/2 

00  30  J-1,LL,K1 

RRIJMRRIJMIRRIJ)  **2“  AL*8E  M*  .500 
30  RRI  J*K’MAL*BE/RR<  Jl 

AO  KMK2 

RETURN 
ENO 


SU5RJUTINE  POISB  I G, $1 , S2, Al , A2, A3, B1 ,B2 ,53, HI , HK, II ,KK I 
C •••  POISSON  EOUATI ON  SOLVER  ••• 

IMPLICIT  REAL*S( A-H.O-Z ) 

C 

REAL  6, SI, S2 

: FIXED  BV  03R.  Pill M-P 1 1 1 1 1 

C 

COMMON  /N5/  RRIAAKLIM 

OMENS  I ON  Gl  1 1 ,1  ,KK!  , SIIII.I.KKI,  S2(II,1,KKI 


0 MENS  I ON  AIIIII,  A2IIII,  A3 1 1 1 1 , BlIKKl,  B2IKKI,  B3CKKI,  MI  IKK)  J 


3*» 


k 


c 


f 


c 


c 


c 


« 


« 


1.  MKt  KK ) 

J 

11 

J-l 

J 

12 

N-l 

J 

13 

Ill-tt-1 

J 

14 

1 12*  1 1*2 

J 

15 

KK1-KK-1 

J 

16 

KK2-KK-2 

J 

IT 

03  10  I «2 . 1 1 1 

J 

18 

G(l,l,ll-GII,l,2l 

J 

19 

GII,l,KKI-G(l,l,KKl) 

J 

20 

1) 

CONTINUE 

J 

21 

03  140  I'liLIN 

J 

22 

00  20  1-2,111 

J 

23 

90  20  N-2.KK1 

J 

24 

S2I  l » 1 ,KI • IRR 1 1 1-831 Kl 1 *GI  I,l,<l*31(KI*G(l,l,K*ll*B2(K)*GI  I.i.K- 

J 

25 

11I-S1II, l.KI 

J 

26 

23 

CONTINUE 

J 

2T 

B991-*4SI2)*RRtl)-42(2) 

J 

26 

34161- 1, 0/B8S1 

J 

29 

4II2I  — 41(21*08681 

J 

30 

00  90  K-2.KK1 

J 

31 

GI2, 1,0-521  2,  l.K)  -04681 

J 

32 

00  SO  I-S.II2 

J 

33 

849**41 1 1 1 *R  R 1 L 1 

J 

34 

0EN4-1. 0/(848* 40111*411 1-1)1 

J 

35 

Mil  I S--41 1 II-0EN4 

J 

36 

33 

GII  .1  .-4S-IS2I  I ,1.KI»42(  I l*G(  1-1.  l.K)  )*OEN» 

J 

3T 

8942- » 1)1  III  )*RR(U*41(  I 11  ) 

J 

3B 

il  111 . l,K)-(  521  I 11.  1,4)  *42(1  1 1I-GI  1 12,  l.K)  ) /(  8932*421  ! 1 1 ) • MI  1 1 12 

J 

39 

111 

J 

40 

S2III l.l.KI-GI III. l.K) 

J 

41 

30  40  IL- 2,  1 12 

J 

42 

1 • 1 12- I L*2 

J 

43 

S 21  I , 1,0 -G<  I,  l.K) -Nil  II-S2I  1*1,  l.KI 

J 

44 

43 

CONTINUE 

J 

45 

SO 

CONTINUE 

J 

46 

00  43  1-2,111 

J 

47 

33  43  K-2.KK1 

J 

48 

Gil, 1,0-52(1. l.K) 

J 

49 

•40 

continue 

J 

50 

33  T 3 1-2,111 

J 

51 

Gl  I , N , 1 1 -Gl 1 , l ,2  1 

J 

52 

Gl  I,N,K<I»G( 1 , 1.KK1) 

J 

53 

TO 

CONTI NJE 

J 

54 

00  80  K-2.KK1 

J 

55 

0(1,1, <1-3(2, l.KI 

J 

56 

GUI, 1,<I—  Gl  I 1 1 , 1 , K ) 

J 

57 

40 

CONTINUE 

J 

58 

33  43  N-2.KK1 

J 

59 

30  93  1-2,111 

J 

60 

S2(I,l,O-(RRIL)-4  3IIII*0l  1.1,  <1  ♦41111-SI  1*1,1,  K)*42  (II  *G(  1-1,1, 

J 

61 

1 K I— SI  1 1 , 1 ,K 1 

J 

62 

TO 

CONTI NJE 

J 

63 

644l«*lil2l*RR(ll-B2(2) 

J 

64 

08381-1.0/8981 

J 

65 

4KI2I — 91(21-08881 

J 

66 

00  120  1-2,111 

J 

67 

Gl  I, 1.2I-S2I  1.  1,21*08881 

J 

68 

00  100  <-3,KK2 

J 

69 

948- *99(KI*R4IL) 

J 

70 

3EN4-1.3/IB49*B2(K)*MKIK-1  II 

J 

71 

WKIKI  — 91IKI*0EN4 

J 

7? 

100 

Gl 1,1,0-1521 1 , 1,KI»B2(KI*GI I, l.K- 11 l*JEN4 

J 

73 

8842- * 9 1 1 KK1 I*RR(L)-B1(KK1 1 

J 

74 

Gl 1 . 1, 5K1I-I  $211 ,1,KK1I*B2IKK1I*GI I ,1,KK2I 1 /( BRB2*R2 I KK l I* WK I KK2 

J 

75 

111 

J 

76 

S2I  I, 1,<K1I-G( I.1.KK1I 

J 

77 

00  110  KL-2.KK2 

J 

78 

K-KK2-<L*2 

J 

79 

52(1,1,0-G(  I,l,KI-4K(KI*S2t  l.l,K*ll 

J 

80 

110 

CONTI  IUE 

J 

81 

120 

CONTINUE 

J 

82 

00  ISO  t-2,111 

J 

83 

>3  113  <-2,KKl 

J 

44 

Gl 1 , 1 ,<I-S2I 1 , l.KI 

J 

85 

ISO 

CONTINUE 

J 

86 

30  140  t-2,111 

J 

87 

Gl  1 , 1 ,1 1 -Gl  1 • 1 ,2  1 

J 

88 

Gil  ,1,<0«GI  1 , 1.KK1I 

J 

49 

140 

CONTINUE 

J 

90 

35 


00  150  K-2.KKI  J 

G(1,N,<CGI2,N,KI  J 

5III,N,<1  — G(  I U.N.KI  J 

150  CONTI  NU5  J 

160  CO-NT  IN'J  5 J 

RETURN  J 

6 NO 

SUB*  JUT  IN 6 TSTEP  I ZU,  ZN  , P,  OXG,  DZG,  10,  JO,  S I GN)  K l 

C •••  NOOtFV  FLOW  BY  PRESSURE  GRADIENT  ••*  K 2 

INPLICir  REAl.«8(  A-M,  0-2  I K 3 

C K A 

REAL  Z'J.ZW.P  K 5 

C K 6 

DIMENSION  ZUIIO.JOI.  ZNllO.JO),  PUD, JO),  OXG(ID),  OZGIJO)  K 7 

C KB 

II- ID  K 0 

JJ-JO  K 10 

III- 11*1  K ll 

112-11*2  K 12 

JJ1-JJ*1  K 13 

JJ2-JJ-2  K 14 

C K 15 

C K 16 

00  10  J-2.JJI  K 17 

00  13  1-2.1! 1 K 18 

ZUIl,  il -ZUII .J)-DXGt D-IPI Cl.Jl-PI I,J)I-S1GN  K 14 

10  C3NTINJE  K 20 

C K 21 

00  20  J-2.JJ2  K 22 

00  20  1-2. Ill  K 23 

ZNI  I,  JCZWII  ,J)-OZG<  JCI®I  I,J*l)-PII,jn*S!GN  K 24 

20  CJNTINJ*  K 25 

C K 26 

00  30  1-2,112  K 27 

ZUII, ll-ZUI 1,21  K 28 

ZUI  I,J  JI-ZUI  I.JJll  K 29 

JO  CONTINUE  K 30 

C K 31 

C K 32 

DO  60  J-2.JJ2  K 33 

C K 

IP  IZUIUJI)  50,50,43  K 

C K 

43  CONT IN  JE  K 

C K 

ZWI1.JC0.0  K 

ZNI  II  , JCZMI  I I 1,  J)  K 

C K 

50  CONTINUE  k 

C K 

ZNl  1,  Jl  -ZWC2  ,J  » K 

ZRI  tt,JI-ZM(  IIl.J)  K 

C K 

60  CONTINUE  K 

C K 

C K 

RETURN  K 

ENO 

SUBROUTINE  BOUNO  IU,H,Q,T«ZBL, ZBR ,ZWL«Z4R,B,BF,It,JJ,DZl)  L 

C •—  APPLIES  BOUNOARY  CONOITIONS  TO  U,J,Q  AND  T •*•  L 

IMPLICIT  REAL* 81 A-H, 0-Z 1 L 

CONN  JN  /N2/  COF.RLAM.PI  ,GA  P , DE  P,  CSTRAT  L 

DIMENSION  UIII.JJ),  WIII.JJ),  TIII.JJI,  3III.JJ),  ZBLIJJI,  ZHLtJ  L 
1JI,  ZSRIJJ),  ZJRIJJI,  BIJJ),  BFIJJI,  0Z1IJJI  L 

REAL  U.4,a,T  L 

JJl-JJ-l  L 

111- 11*1  L 

112- 11*2  L 

ZERO- 3 L 

Jill  1,1  CUC I 1 1, 2 I L 

Ul  1 II  .JJI-UI  II  1,  JJ1)  L 

Til  11,11— Tl 1(1,21  L 

TIIll.JJI-TIIll, JJ10  L 

C L 

C SIDES  L 

l 


DO  10  J-2.JJ1 


on*-  no  o n 


« 


L 

LEFT  L 

IF  IZDLIJI.Ej.O.I  QI1.JI*)I2iJI  L 

IF  IIIU  JI.EQ.O.)  T(l,JI-T(2,JI  L 

IF  (ZMLl JI.EQ.0.1  W(  l.J)*M(2,JI  L 

IF  (ZBLUI.LT.O.I  T( l.JI«T(2.JI*2.-T(3,J)  L 

IF  (ZBLUI.LT.O.I  Q ( l , J I *Q( 2 , J 1*2.-013, J I L 

IF  (2SLUI.LT.0.)  Q<l,JI«ANAXl(0(  1 , J l ,Q( 2, J I /4. 1 L 

IF  (2KLUI.LT.  0.  I W(l'J)«W(2tJ)*2.-W(3>J!  L 

L 

RIGHT  L 

a-dmax i ( ( ( u(  1 1 ii  j* 1 1 in,j-in/32i(j»/2i**2-(Tun, j»n-Tiin.  l 
1J-UI/0ZIU)/2,ZER0>  L 

IF  ( 2 3R ( JI . GT .0.1  01 II • JI-BF( JI-A-RLAM-RLAM/COF  L 

IF  (ZSRU).GT.O.)  T(II,JI-a(JI  L 

IF  (ZM1I JI.5T.0. I W(II,JI-0.  L 

IF  (ZBKJI .GE.O.I  U< II, Jl-U( 112, JI  L 

IF  (2MRUI.EQ.0.I  m(II,JI-W( Ill.JI  L 

IF  (Z3R(  JI.EO.O.I  Qdl,  Jl-Odll,  JI  L 

IF  (ZBRU  I.EQ.O.I  T( II, j|>Tllll, JI  L 

IF  (Z3R( JI.LT.O.)  T( II , J)«T( III, JI-2.-TI 1 12, JI  L 

IF  (2BRUI.LT. 0.1  0(11.  JI-QI  III,  JI-2.-0I 1 1 2, JI  L 

IF  I23RUI .LT.O.I  0( II, J)*AM4X1(0( II, J t.OI I ll.JI/4.1  L 

IF  (2  81 ( J I «LT. 0. I Ullli  JI-JI  Ml.  J I *2  .-J  (I  I 2 . J I l 

IF  (ZM1( JI .LT .0.  I M( II ,JI-R( II l, J)*2.-M( 112, JI  L 

0 continue  l 

L 

TOP  AND  80TTH  L 

00  20  t-l,  II  L 

0(1  .ll*J(  1,21  L 

T(I,1I-TI  1,21  L 

J ( I . 1 I * J ! 1 .2  I L 

J(l.JJI*U( I.JJ1I  L 

r ( 1 1 j j i • r ( i.jjii  l 

0(1 ,J JI »0( I, JJ1I  L 

20  CONTINUE  L 

RETORN  L 

ENO 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


t 


SUBROUTINE  OUTPUT  ( U.  2 J . X.  2W  ,T  , 0 , F , E , $ , 40V , R , X , 2 ,Y , DX  1 , 021  . RL  41,  M 
10T.TINE, J4P.0EP,  II .JJ.NC, NR.NXN, 12 . J 1 ,8, 0X2. 022. TT ,TB, 23.0SCALE.NN  M 


2.TIN,M10FT,G,H,NC0PYSI  M 

C OUTPUT  ROUTINE  CALLED  EVERY  TINE  STEP  •**  M 

IMPLICIT  R E AL- 8(  A-H,  3-2  I M 

COMMON  /HZ/  COF,XXXX.PI,VYY,222,CSTRAT  M 

COMMON  /TXTC/  I T XT ( 60 1 . 0 AY . J TX T ( 60 1 M 

COMMON  /STREAM/  IS  M 

REAL  U.ZU.W.ZW.S.P.ZQ  M 

REAL  T, J.E.F .AOV.G.H  M 

DIMENSION  Fdl.JJI,  EdI.JJI,  Odt.JJ).  Mdl.JJI  M 

DIMENSION  U(II.JJ).  Rdl.JJI.  T ( 1 1 , J J I . Odt.JJ).  Sdl.JJI,  Rltl  M 
II.  Xdll.  2 ( J J I . YIJJI,  OXUIII,  0 2 1 ( J J I . ZUdl.JJi.  ZMdI.JJ).  3(  M 
2 JJI.  0X2(  1 11.  02 2( JJ I . 20(11. JJ)  M 

DIMENSION  AL ( 4 1 * AR(4I,  ICMI.I,  AOVIJJ.III  M 

OAT  4 ICM/1HT.1HU.IHA.UT/  M 

0 1 ME  NS  I ON  F l ( 7 5 1 . F2( 75  I , e3(7il,  F4(75l,  F5(75l,  F6( 751  M 

0 1 MEN  SION  IFR0M(3I,  ICNF  R(  31  M 

DIMENSION  SX ( 8 I . T IN ( 2 1 , TXI3I  M 

DIMENSION  I T SAR ( 14  I , TAR(SI  M 

DATA  I T JAR  / l, 2. 3.4. 6. 8,  12.  18,2  4.  36.48.72.36.120/  M 

C M 

Itl-ll-1  M 

II2-1I-2  M 

JJ1-JJ-1  N 

JJ2-JJ-2  M 

T0-TT-T3  M 

C M 

00  10  1*1.8  M 

10  SXIII-3.  M 

SX(6l-l.0l0  M 

S4( 71*1.010  M 

S X ( 81*1.010  M 

SUM7-0.J  M 

C N 

00  40  J-1.JJ1  M 

DO  40  1-1. Ill  M 

Txm-;i  i.ji  m 

TXI2I  ■ )(  I , JI  M 

TXI3I-2K  I.JI  M 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 
2b 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


i 


c 


37 


i 


t 


( 


€ 


t 


» 


ft 


ft 


C 


ft 


ft 


ft 


TX(  4 1 »'J(  I f J) 

M 

41 

TX( 51-4(1. J) 

M 

42 

TX(4)-3A8S(TX(4I  1 

M 

43 

T X( 51 -3AdS(TX( 51  ) 

M 

44 

03  20  (-1.5 

M 

45 

20 

SXI X) * 3 MAX 1(SX(K),TX(KI) 

H 

46 

30  50  <>6.8 

H 

47 

SXKI- JNINH SX(K»,TX(X-5)) 

M 

48 

53 

CONTINUE 

M 

49 

*3 

CONTINUE 

M 

50 

03  SO  J-2.JJ1 

M 

51 

30  S3  1-2,111 

M 

52 

o-( 3( i , j i -ui  i-i.j) i/oxi iiaiwa 

,j)-w(i,j-i)i/ozkj) 

M 

53 

SJM7-S JN740-0 

N 

54 

53 

CONTINUE 

M 

55 

0 1 V-SlMT 

M 

56 

C 

M 

57 

C 

M 

58 

T INI  11-3 

M 

59 

TINI21-3 

M 

60 

A -0 

M 

61 

C-0 

M 

62 

IT-1 

M 

63 

JT-O 

M 

64 

OJ  70  J-2.JJ1 

M 

65 

Ie  (U( l.JI.CE .0. 1 GO  TO  60 

M 

66 

JT-l 

M 

67 

A-A-OI  1.J1-I T( 2, Jl -1.5-T (3,J  )-. 

51-OZK  Jl 

M 

68 

C-C  *U( 1.JI-32U Jl 

N 

69 

IF  (J  .S3.JJ1  1 GO  TO  60 

M 

70 

GO  TO  73 

M 

71 

63 

IF  (JT.cQ.O)  GO  TO  73 

M 

72 

TIN(IT1-A/C*T0*TB 

M 

73 

JT-0 

M 

74 

A-3 

M 

75 

C-0 

M 

76 

I T»  I T*1 

N 

77 

IF  I IT  ,GE  .5  • GO  TO  80 

M 

78 

73 

CONTINJc 

M 

79 

33 

CONTINUE 

M 

80 

C 

M 

81 

IF  (NK.EQ.2I  PRINT  90 

M 

82 

IF  (NC.SQ.O)  PRINT  90 

M 

83 

JO 

FORMAT  1/*  NC  TMAX  3 MA  X ZQMAX  UMAX 

WMA  M 

84 

IX  T MI N QMIN  2QMIN  OIV  TIN1 

TIN2  M 

85 

2 */ 1 

M 

86 

c 

M 

87 

PRINT  130,  NC.SX.OIV.TIN 

M 

88 

100 

FORMAT  ( 16, IP, 11011.4) 

M 

89 

IF  (MK.SQ.NKMI  PRINT  93 

M 

90 

2 

M 

91 

c 

M 

92 

IF  (NK-NKMI  480,110,113 

N 

93 

c 

M 

94 

c 

M 

95 

ao 

C3NTINJE 

M 

96 

NK-0 

M 

97 

c 

M 

98 

c 

M 

99 

C WRITE  ANSWER!  AT  ENO  OF  STREAM  13  F3R 

PRINTING 

M 

100 

C 

M 

101 

REWIND  10 

M 

102 

READ  113,120) 

M 

103 

120 

F3RMAT  ( 29(  / ) ) 

M 

104 

C 

M 

105 

IF  (TIN(l).EQ.O)  GO  T3  150 

M 

106 

IF  (T  I((2).EQ.0t  WRITE  (13,1331 

TfN( 11 

M 

107 

IF  1 T I N(  2 ) .NE. 0)  WRITE  (13,143) 

TIN 

M 

108 

150 

FORMAT  IIOX.'MEAN  PLANT  INFLOW 

TEMPERATJRE  IS'.FS.l) 

M 

109 

140 

FORMAT  (6X, 'MEAN  PLANT  INFLOW  TEMPERATURES  ARE*,2F6.1I 

M 

110 

150 

CONTI  N JE 

M 

111 

C 

M 

112 

c 

M 

113 

C SCALES  IN  FEET 

M 

114 

00  163  1-1,11 

N 

115 

XII  >-<( I l-OSCALE 

M 

116 

R( I I* R ( I l-OSCALE 

M 

117 

163 

CONTINJE 

M 

118 

00  170  i-l,JJ 

M 

119 

IUI-21  J l-OSC  ALE-OSC  ALE 

M 

120 

1 1 


c 


r 


r 


r 


c 


c 


c 


c 


« 


c 


VIJI-YI JI*DSCALE-OSCALE 

H 

121 

170 

CONTINJE 

M 

122 

C 

H 

123 

C XPLOTS 

H 

12* 

C 

H 

125 

PRINT  ifO 

M 

126 

ISEP-1 

H 

127 

JLINE-Z5 

M 

128 

C SNVl  RONNENT  AL  IMPACT  PLOT 

H 

129 

JSC-0  ilRT(DSCALE*T0-32*l.3E-*) 

M 

130 

F2I1I-J 

M 

131 

00  180  J“i,JJl 

M 

132 

F1(J|»IT(  II,  J)*T(II,J*1)  1/  2.  *T0*TB 

N 

133 

180 

F2(J*ll«F2(J)*U( II l»J+l)*0Z2tJ*l)*DSCALE*4I0FT*USC 

M 

13* 

TRT»F 1 C JJ1  ) 

N 

135 

TR9-FII11 

M 

136 

CALL  X»L0T  (F2,F1,JJI,1,J0,TR3,TRT,'CUFT/S  ADDED  BELOW  T **','T 

M 

137 

1 INE  «•  ,T  ME  1 

M 

13R 

C 

M 

139 

C 

M 

1*0 

XL-0. 0 

M 

1*1 

XJ- GAP*  JSC  ALE 

M 

1*2 

C 

M 

1*3 

DO  190  1-1,1 I 

M 

l** 

Fill (»T( I , JJ1) *TD»TB 

M 

1*5 

190 

CONTINJE 

H 

1*6 

CALL  X BL0T  IF1  ,R  ,1 I , ISEP, JLINE.XL.XU, • T ( I , JJt)  ', 

H 

1*7 

I'TME  *'  .TIME! 

M 

1*8 

c 

H 

1*9 

C DOUBLE  PRINT  COUNT  AT  END. 

M 

150 

!P*l 

H 

151 

IF  (NC.GE.NN)  IP-2 

H 

152 

00  250  IPC-l.IP 

M 

153 

IS-6 

M 

15* 

PRINT  A 70 

'I 

155 

NCONT-12 

M 

156 

IF  ( IPC.E3.2  1 NCONT-20 

H 

157 

ISEP-l 

M 

158 

NLIN--1 

* 

159 

4NAG-1 

M 

160 

2ER- 0 

M 

161 

AU-GA®*0SCALE/ IPC 

M 

162 

ZL--0SJALE 

M 

163 

c 

M 

16* 

CALL  CCNTOR  ( 0 ,R  ,2 , 1 I , J J , NCONT , I SE  P,  NL I N , A NAG.  Z ~ P . XU.  ZL  . ZE  R , • TUR 

M 

165 

1BULENCE  ISOLINES  ••  ’ « 1 T 1NE  ’.TIME) 

166 

c 

M 

167 

CALL  CO  (TOR  (F,R  ,2  , I I , J J .NCONT , I SE  P.  NL  I N , AMAG.  Z ER . XU,  ZL , ZE  R , • TUR 

M 

168 

1BULENT  OlFFUSIVITY  ••'.•TINE  *,TINE) 

M 

169 

c 

M 

170 

c 

M 

171 

DO  200  1-2,111 

M 

172 

00  200  J-2.JJ1 

M 

173 

200 

Ell, J)«CSTRAT*RLAN*RIAM*(TII,J*1I-T(I,J-1) )/2./DZl(JI/Q(  t,J) 

•4 

17* 

CALL  10  ( E , 1 I. JJ) 

H 

175 

c 

- 

M 

176 

CALL  CO  (TOR  (E  ,R,Z,I!,JJ,NCONT,ISEP,NLIN,AMAG.ZFR,XU,Zl,ZER,'CST 

M 

177 

1 R AT  • N**2  • L *• 2 / 9 '.'TIME  '.TINE) 

M 

178 

DO  220  1-2,1  11 

H 

179 

DO  210  J-2.JJ1 

M 

180 

AA«(U(  I.JI-JI  I-1,J)  l/OXKI  1 

N 

181 

A8-1  J(  I , J*1)*JII-1 , J»1)-JI  l,  J-l)-Jt  1-1,  J-l  ) )/DZl  1 J 1 ♦!  Ml  1*1  ,J)*W( 

M 

182 

11*1,  J-l  1-41  I - 1 , J 1- W(  1-1,  J-l)  l/OXK  II 

M 

183 

SI  I ,J>**.*AA*AA*AS*A8/16»1.0-10 

N 

18* 

A-SI  I , J 1 

M 

185 

El  I,  JI-OSORT ( A)/ (FI  I , J ) ♦l.E-ZOI 

M 

186 

SI  I,  J)-(  T(  l,J*l l-TI I ,J-1)) /2./0Z1IJI/SI I ,J  1 

* 

187 

210 

CONTINJE 

H 

188 

220 

CONTINJE 

M 

189 

C 

190 

CALL  BO  (S.II.JJI 

M 

191 

CALL  BO  1 E , 1 I • JJ ) 

N 

192 

c 

* 

193 

CALL  CC  ITOR  IS  ,R  ,Z,  I I.JJ, NCONT, I SE P , NL 1 N , AN AG, ZE R, XU, ZL , ZER, 'RIC 

N 

19* 

1NAR0S0N  NUMBER  ••  '.'TIME  ',TIMEI 

* 

195 

c 

M 

196 

CALL  CCNTOR  IE.R.Z, 11, JJ,. NCONT, ISEP, NLIN, AMAG. ZER, XU, ZL, ZER,'  RE 

* 

197 

1VNOLOS  < JNBER  ••  ','TIME  '.TIME! 

H 

198 

c 

H 

199 

00  2)0  J-l.JJ 

H 

200 

39 


1 


c 


€ 


€ 


€ 


t 


C 


C 


C 


C 


C 


C 


30  233  1-1,11 

M 

201 

233 

Stl,J)»TB«-TO*T  ( l,J  1 

N 

202 

CALL  COITOR  ( S ,R  ,1 , 1 1 , J j ,NCO  NT  , I S E P,  Nil  N,  AM  AO,  2E  R.  XU,  2L  , ZER  , • TEH 

M 

203 

IPERATURc  •*  •.'TIME 

•.TIME) 

N 

20* 

C 

M 

205 

c 

M 

206 

00  2*0  1-1,11 

M 

207 

Eli, 1)0.0 

K 

208 

SUN-3.  3 

M 

209 

00  2*3  J-2.JJ 

M 

210 

SUM-SJN»UI  I , J)  -DZ1  ( J) 

M 

211 

El  I, J l-SUH 

M 

212 

2*0 

CONTINUE 

N 

213 

c 

M 

21* 

CALL  CCNTOR  1 E . X ,T , 1 1 , J J , ICO  IT  , 1 SEP , NL I N , A MAC, ZE « , XU , 2L , ZE R , • STR 

M 

215 

1 E AH  LINES  •*  • » • T IMe 

• .TIME) 

N 

216 

C 

N 

217 

C 

M 

218 

250 

CONTINUE 

M 

219 

C 

m 

220 

C 

M 

221 

C VPL3T 

RANGES. 

M 

222 

I-SX(*)*12*USC+.9 

M 

223 

uscp»i/io. 

M 

22* 

jscM»-uscp 

M 

225 

00  263  J-1,JJ 

H 

226 

2 50 

P1(J)-T(1,J| 

N 

22  T 

CALL  N A <N IN  ( F l , Z,  J J , TMX  ,T  NN  , 

TAB ,ZMX .ZHN.ZAB) 

N 

228 

IF  ITMX.LT.l)  TMX-1 

M 

229 

IF  ITNN.GT.3I  TMN-0 

M 

230 

ITL«TB*-T0-TMN*.01 

M 

231 

I T'J«T3*TD*TMX*.99 

M 

232 

ITO-I TJ-ITL 

M 

23? 

30  270  1-1,1* 

M 

23* 

I T S-I  T 5 ARC  I ) 

N 

235 

IF  < IT3.LE .ITS ) GO  TO  290 

M 

236 

270 

CONTI N JE 

H 

237 

2 30 

CONTINUE 

n 

238 

1TL-ITL-I  I TS-  I TO  )/  2 

M 

239 

ITU-ITLMTS 

M 

2*0 

DO  293  1-1,5 

M 

2*1 

290 

TARll )*1TL*1TS»1 1-1)/*. 

M 

2*2 

z 

N 

2*3 

C VPLOT 

CALL  LOOP. 

M 

2** 

00  3*3  1-1,111,112 

M 

2*5 

30  300  J-l.JJ 

N 

2*6 

A3V(J,1)-T(I,J) 

M 

2*7 

AJVI J,2)»J( I ,J) 

N 

2*8 

A)V<J,i)«B(  J) 

M 

2*9 

30  J 

CONTINUE 

N 

250 

C 

N 

251 

ALI 2) --NSCP/USC 

M 

252 

AR  ( 2 ) -UiCP/USC 

M 

253 

ALI  31  -It  TL-T3)  /TO 

M 

25* 

ARI 3) *1 ITU-TBI/TO 

M 

255 

ALI l)- ALI 3) 

M 

256 

ARI  11  -ARI 3) 

M 

257 

c 

N 

258 

00  333  lPC-1,2 

N 

259 

IS- 9 

M 

260 

IF  ll.EJ.ll  IS-8 

M 

261 

IF  IIPO. E3. 21  IS-IS-3 

M 

262 

CALL  7 ’L OT  1 AOV,  JJ  , 3 ,2 , Z L,  ZER 

, ICH,AL,AR,5*!PC,7,*-IPC,6, IS ) 

M 

26  3 

IF  1 1 9C.E0.1 ) WRITE  IIS,  113) 

TAR.USCM.USCP 

M 

26* 

IF  IIPO. EO. 2)  WRITE  IIS, 323) 

TAR.USCM.USCP 

H 

265 

310 

FOR* AT  1/'  0EG',F6.1,3F6.1,F5 

. 1/ • FT/S', F5. 1,9X, *0. 0* ,F11.1) 

M 

266 

320 

FORMAT  l/«  0EG',F6.1,3F12.1,F11.1/'  FT/ S • , F5. 1, 2 IX , • 0. O' , c 23. 1 ) 

N 

26  7 

33  3 

CONTINUE 

M 

268 

3*3 

CONTI  NJc 

M 

269 

C 

M 

270 

1FR0NIII-3 

N 

271 

IFRONI >1-10 

M 

272 

IFROMI 3 ) -9 

M 

273 

ICHFRI 11-36 

M 

27* 

ICHFR  I i )-60 

M 

275 

ICHFAID-32 

M 

276 

C 

M 

277 

c 

N 

278 

C START 

MULTI *L2  COPY  FORMAL  OJT»JT. 

M 

279 

I S-6 

M 

280 

Ln 


IF  (NC.EQ.NN.AN0.MC0PYS.3r.il  (S>13  H 

C M 

HRITE  I IS, 470)  M 

C M 

NCONT  » ) N 

XU-GAP*  JSCALE  N 

C N 

00  350  J-UJJ  N 

00  ISO  1-1,11  N 

F(I,J)-F(l,JI*USC*OSCALE  m 

350  S(1,J)*TB+TD*T(I,J)  M 

CALL  CQNTOR  (S  »R  ,2 . 1 1 . J J ,NCO  IT  , I SEP ,NLI N , ANAG,  2E R,  XU,  21 , 2ER , • TEN  H 
l PERATUR  E ••  ','TIME  '.TINE!  M 

C X 

C n 

00  360  (-1,11  M 

E(  I,  1 1- J.  0 M 

SUM-0.)  N 

00  160  U-2 , J J N 

SUN-SU4*U(  I.JI-OZK  JI*USC*0SCALE*U1DFT  N 

EII.JMSUM  M 

343  CONTI  Xu!  M 

C M 

Jl-JJl  N 

CALL  C3NT0R  C E ,X  ,Y , 1 1 . J 1 .NCONT  , I SEP , NLI  N,  ANAG,  IE  R,  XU , 2L,  2ER , • STR  M 
IEAN  LINES  (CUFT/SI  •••,'TIHE  • , T I HE ) M 

C H 

CALL  C3H81N  (42, IS, 132, 3,1 PROH, ICHFR I M 

C END  FIRST  P 43c  PAIR  M 

C H 

HRITE  I IS, 4701  M 

C N 

CALL  CONTOR  (F ,R ,2 , 1 1 , J J , NCONT , I SEP , NLI N.AMAG, IE R, XU, 2L, IER , • OIF  H 
IFUSIVITY  (FT**2/SI  •*• , ' T ME  • , T ME)  N 

C N 

HRITE  ( IS, 3701  M 

370  FORMAT  1111,'M  T E MPE  A ATURE  T 4N0  H0RI20HTAL  MOTION  U AT  LEFT  4 M 

1  NO  RIGHT  BOUNDARIES,  WITH  THc  ANSI  ENT  TEMPERATURE  A »*•/)  M 

1 FROM  111-11  M 

I FROM (21-12  M 

ICHFRID-72  M 

ICHFR (21-96  M 

CALL  C31BM  (75,  IS, 132, 2, 1 FROM, ICHFR)  M 

C M 

REMIND  10  M 

00  390  1-1,11  N 

READ  I 10,4101  TX  M 

380  HRITE  (3,4101  TX  M 

00  390  1-1,17  N 

READ  (10,410)  TX  M 

390  HRITE  19,410)  TX  M 

00  400  t-1,6  M 

REAO  I 10,410 , END-420  I TX  M 

400  HRITE  I 1,410)  TX  N 

410  FORMAT  I IX, BAS)  M 

420  CONTINUE  N 

C M 

1 FROM!  1 )-S  N 

I FROM  I > ) -9  M 

ICHFR III -68  M 

ICHFR I 21-60  N 

CALL  COMBIN  (20, IS, 132,2, IFR3H, ICHFR)  M 

C M 

IF  (IS.EQ.6I  GO  TO  440  N 

00  430  I-i.NCOPVS  M 

430  CALL  CO  IB  IN  (289,6,128,  1, 1 3, 129 ) M 

1 S*6  M 

440  CONTINUE  M 

C M 

PRINT  470  M 

C END  MULTIPLE  COPY  FORMAL  OUTPUT.  M 

C M 

C RETURN  TO  OMENS  IONL ESS  SCALES  M 

00  490  1-1,11  M 

XIII-XI  D/DSCALE  M 

A ( 1 1- A(  t l/OSCALE  M 

490  CONTINUE  M 

03  440  J-l,JJ  M 

Y(  JI-VI J) /OSC ALE •!  N 

2 ( J) -2(  J 1 70SC  ALE  *1  M 


281 

282 

283 

284 

289 
286 

287 

288 
284 

290 

291 

292 

293 

294 
299 

296 

297 

298 

299 

300 

301 

302 

303 

304 
309 

306 

307 

308 

309 

310 

311 

312 

313 

314 
319 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 
32  8 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 


CONTINUE 


r 


c 


€ 


C 


€ 


( 


460 

C 

470 

C 

440 

C 


FORMAT  UH1////) 

CONTINUE 

RETURN 

END 


N 362 
N 363 
N 364 
N 365 
H 366 
N 367 


SUBROUTINE  CONTOR  ( Q,XV , Y9, NX, NY .NCONT , I SEP ,NL IN , AMAG,XL,XU, YL, V N 
1U,  TITLE,  iJBTIT.VARl)  N 

C N 

C GL  Y Y*  S LINE  CONTOUR  ROUTINE  ***••  N 

C 0(1, J)  IS  VALJE  AT  (XV (II, Y3 (Jilt  WITH  XV  ANO  TO  MONOTONIC  INCREASING  N 
C DIMENSIONS  ARE  NX  ANO  NY  N 

C APPROXIMATELY  NCONT  CONTOURS  ARE  DRAWN  N 

C OUTPUT  SPLIT  INTO  ISEP  PARALLEL  PAGES  N 

C VERTICAL  SCALE  DETERMINED  BY  APPROXIMATELY  NL IN  HORIZONTAL  ROWSt  N 

C IF  NLIN.LT.O  , PHYSICAL  SCALINGt  WITH  VERTICAL  STRETCHED  BY  AMAS  N 

C CONTOURS  REGION  XL. LE .X .LE.XU  , YL .L E. Y.LE .YU  N 

C TITLE  IS  24  CHARACTERS  AT  8 PER  WORD.  N 

c subtit  is  8 characters,  e.g.  *time  - • , with  vari  the  time  value.  n 

C MAIN  PROGRAM  SHJULO  CONTAIN  — N 

C COMMON  /rxTC/ITXT(60l,0AY  N 

C RE  ADI  5,  1 1 ITXT  N 

Cl  FOAMAT(60A1I  N 

C CALL  I)AY(OAY)  N 

C --  TO  READ  FJRTHER  TEXT  OFF  A CARO  AMO  SET  OATE  IN  OAY.  N 

C BLANK  COMMON  IS  USED  FOR  WORKING  SPACE  1 962  WORDS!  AND  SO  PASSED  N 

C ARRAYS  SHO'JL)  NOT  8E  IN  BLANK  COM  ON.  N 

C N 

IMPLICIT  REAL*8(A-HtO-ZI  N 

C N 

DIMENSION  IC(120.2I  N 

DIMENSION  LGII21I.  CHAR (601,  TITLEI3I,  L0I121I  N 

CUMMON  /TXTC/  IT  XT (601,0 AY  N 

COMMON  (CtLO.LC.G, H N 

COMMON  /STREAM/  IS  N 

INTEGER  BLANK, MINUS, PLUS, BAR, CHAR  N 

DIMENSION  XV I NX)  , YOINY),  Gl  120,4),  XK11I  N 

REAL  F,  ,,  I C,  H ( 4 1 N 

REAL  0( MX, NY)  N 

YVI J I -YD ( J I N 

F(  l ,JI  . )(  I ,J)  N 

C STATEMENT  F'JYCTION  CAN  8E  USED  TO  TURN  A OVER  N 

DATA  BLANK/IH  / , MI NUS/1 H-/ , PLUS/ 1H»/  N 

DATA  BAR/1 H| / N 

C OATA  9AR/1HI/  N 

DATA  CHAR/IHT,  1HS,  l MR,  1 HD,  LHP,  1H0,  IHN,  1HM,  IHL,  1HK,  IHJ,1HI,  IHH,1H  N 
1G, 1HF,1H=,1H0, IHC, 1H8,1HA,1H0,1H1, 1H2, 1H3, 1H4, 1H5, 1H6, 1H7 , 1H8 , 1H9,  N 
21HD,  1H1,  1H2,  1H),  1H4,IH5,1H6,  1H7,  1H 3,  1H9, 1H0 , 1H1 , 1H2, 1H3 , 1H4,  1H5 , 1H  N 
36.1H7.1H3, 1H9, 1H0, 1H1.1H2 , 1H3, 1H4, 1H5, 1H6, 1H7, 1H8, 1H9/  N 

YG-YU-YL  N 

C •••**•*•*•*»«•*••*•••*•••*•••*••*•**•••••••••••*•*»»•*•*•*•*•*****••*«  N 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


IF  (ISeP.GT.4)  ISEP«(!SEP*l)/4*4 
C 

XG*(XU-XL)/ISEP 
00  390  ISP-1,1  SEP 
XB-XL ♦<G*I SP-XG 
00  10  1-1,121 
L0( I )-)LANK 
10  LG(  I l-MI  NUS 

DO  20  t -1 , 121 , 12 
20  LGI 1 1 -PLUS 

00  ID  (-1,11 

30  X1III-A4»XG*II-1)/10 

C 

C FIND  MAXIMUM,  MINIMUM,  AND  CONTOUR  ORIGIN  ANO  INCREMENT 
FMX-F (1,1)  < 

FMN-FMX 
00  40  J-l.NY 
DO  40  I > 1 , NX 
FV-F(I,J| 

FMX-DMAXKFMX.FVI 
43  FMN-OMINl ( FMN, FVI 

IF  IFMX.EO.FMNI  WRITE  I IS, 4001  TITLE, SUBTI T, VAR  I ,FNX 
IF  I FMX. EQ.FMN)  RETURN 

C NOTE  FORMAT  CHANGED  FROM  0 TO  G IN  INOEPENOENT  CHANGE. 
A«OLOG10(  I FMX-FMNI /NCONT  >U1S35 


N 46 
N 47 
N 48 
N 49 
N 50 
N 51 
N 52 
N 53 
N 54 
N 55 
N 56 
N 57 
N 58 
N 59 
N 60 
N 61 
N 62 
N 63 
N 64 
N 65 
N 66 
N 67 
N 68 
N 69 
N 70 


I 


I -A 

N 

71 

IP  IA.LT.O)  l-I-l 

N 

72 

B-I 

N 

73 

C-A-I 

N 

74 

FINC-13.P*B 

N 

75 

IF  1C. LT.. 15051  CO  TO  50 

N 

76 

F INC- 10. --B* 1. 5 

N 

77 

IF  (C.LT..3010)  GO  TO  SO 

N 

78 

FINC-13.P«B«2 

N 

79 

IF  (C.LT..4S16)  GO  TO  SO 

N 

80 

FINC-1Q.«*B*3 

N 

81 

IF  (C.LT.. 65051  GO  TO  SO 

N 

82 

FINC«13.**B*5 

N 

83 

IF  1C. LT. .86901  GO  TO  SO 

N 

84 

F INC* 13. ••B*T. 5 

N 

85 

CONTINUE 

N 

86 

J-FMN/FINC 

N 

87 

K-FNX/FINC 

N 

88 

I-J 

N 

89 

IF  IK-J.GT.25)  I-IJ-X>/2 

N 

90 

IF  IK.LT.26.AN0. J.GT.-26I  t-0 

N 

91 

FOR- 1 *F  INC 

N 

92 

N 

93 

IF  (XT. GT. XV (NX). OR.XB.LT. XV (11)  GO  TO 

113 

N 

94 

IF  (VU.GT.VV(NV) .3R.YL.LT. VV(l >)  GO  TO 

110 

N 

95 

Jd-l 

N 

96 

IB-1 

N 

97 

00  60  4-1. NY 

N 

98 

IF  (VV(J).LT.YL)  JB-J-l 

N 

99 

IF  (YV(J).GE.VU)  GO  TO  70 

N 

100 

CONTINUE 

N 

101 

JT-J-1 

N 

102 

XT-XD-XG 

N 

103 

DO  80  I - 1 1 NX 

N 

104 

IF  (XV(I).LE.XB)  IB-l-1 

N 

105 

IF  (XV(I).GE.XTI  SO  TO  93 

N 

106 

CONTINUE 

N 

107 

IT-I-1 

N 

108 

00  100  J-J8.JT 

N 

109 

00  100  I-IB.1T 

N 

110 

IF  I 0A3SIFI J » J )-FOA)/F 1 NC.GE.. 5)  GO  TO 

110 

N 

111 

CONTINUE 

N 

112 

VRITE  (IS. 410)  TITLE, SU»TIT,VARI,XB, XT, 

YLtYU 

N 

113 

GO  TO  190 

N 

114 

CONTINUE 

N 

115 

- 

N 

116 

NL1NE-1LIN 

N 

117 

IL-NL I NE  Z7.-.5 

N 

118 

IF  INLINE. LE. 3)  IL-72-YG/X3-AN43/7-.5 

N 

119 

IL-NAX3I IL, 1 1 

N 

120 

NLI NE -7*1 L 

N 

121 

CALL  SH0RT7  (YU.FR.IX) 

N 

122 

VRITE  (IS, 420)  TITLE,  1 TX T,  SU3T IT , V AR I , FNX, FHN, FOR, F INC.O  AY , XI 

N 

123 

MRITE  (IS, 4501  YU, LG 

N 

124 

N 

125 

SET  UP  G 

N 

126 

00  120  1-1,120 

N 

127 

IC( 1,21-1 

N 

128 

CONTINUE 

N 

129 

NYU- 3 

N 

130 

NXL-1 

N 

131 

00  100  J-l.NY 

N 

132 

IF  IYU.GT.VVIJ))  NYU- J* 1 

N 

133 

00  140  I-l.NX 

N 

134 

IF  IXB.GT.XVI 1)1  NXL-I 

N 

135 

J-0 

N 

136 

JY-NI  N 3 ( NY* 1 »N YU-2 ) 

N 

137 

N4-1 

N 

138 

00  210  N-NM.4 

N 

139 

N-JY-N 

N 

140 

IX-NAX3I0.NXL-2) 

N 

141 

IX-NINOI IX, NX— 4) 

N 

142 

00  200  1-1,120 

N 

143 

X-XB+I !-.5)*XG/120 

N 

144 

IF  (X.LE.XVI IX-3I.0R.IX-4.E0.NX)  GO  TO 

170 

N 

145 

IX-IX-1 

N 

146 

GO  TO  160 

N 

147 

e-o 

N 

148 

00  1«0  K-1,4 

N 

149 

0-F ( IX*R,N) 

N 

150 

03  130  L-l,* 

f 

180 

IF  IL.SO.K)  GO  TO  180 

O-O-IX-XVI IX*L ) ) /( XVIIX*K)-XVI IX*L) ) 
CONTINUE 

190 

2-6*0 

200 

Gt I, Rl-E 

JJ-J 

IF  IJ.NE.O)  GO  TO  220 

210 

C0NTINJ6 

c 

C 

c 

MAIN  PROGRAM 

JJ-1 

220 

00  360  J> JJ, NL INE 

IERC-0 

Y-YU-Y  G-I J-.5I/NLINE 

IF  (T.GE.TVI JY-3I.OR.JT.E3.5)  30  TO 
JT-JT-1 

c 

230 

00  230  K-1,3 

00  230  1-1,120 

Gl I ,K ) *G( I ,K* 1 ) 

MM-* 

GO  TO  150 

2*0 

00  270  1-1,120 

6-0 

00  260  K-l,* 

c 

250 

0-31 I ,K) 

00  250  L-l,* 

IF  (L.68.KI  GO  TO  250 

0-0-1  V-VVC  JT-L  ) ) /I  YVI JV-K)  -YVI  JT-L  ) ) 
CONTINUE 

260 

6-6*0 

ICI 1,11— IC(I,2) 

ICC  1 , 2 ) - 1 E-FOR  )/FINC*22 

€ 

• 

270 

CONTINUE 

IF  IJ.5Q.il  GO  TO  360 

IQ-0 

00  310  1-2,120 

LSI  II-3LANK 

HID-  ICI  1-1,1  I 

HI  21- ICI 1-1,2) 

HI3I-ICI  1,21 

« 

HI  6 1 - I C 1 1,11 

00  280  K>1,3 

N-1*K 

00  230  L-M,* 

IF  IHIU.LE.HI  LI  I GO  TO  28  0 

HH-NIK) 

HIKI-HIL) 

HILI-HH 

c 

230 

CONTINUE 

IL-IHI l )*HI 2)1/2. 

IU-IHt3l*Ht*ll/2.-l 

IF  IIU.LE.60.AN0.IL.GE.il  GOTO  300 
IF  II53C.LT.1)  WRITE  IIS, 2901  IU,IL,I 
IERC-1 

290 

FORMAT  | • IU,I L, I,E-  ',311*. 01*.*) 

GO  TO  310 

« 

300 

CONTINUE 

IF  ML.ST.IU)  GO  TO  310 
LGItl-CHARIILI 

IF  I IL.EQ. IU)  GO  TO  310 

loi 1 1 -Chari iu) 

18-1 

310 

CONTINUE 

C 

c 

C 

PRINT  LINS,  POSSIBLY  WITH  Y VALUE 

i 

JLIN6-J-1 

IF  1 JL  INEZ 7-7.  NE  .JL  INE)  SO  TO  320 
Y-VU-YG-JL INE/NLINE 

CALL  SH3RT7  IY,FR,IXI 

WRITE  IIS, *50)  V ,LG 

GO  TO  330 

1 - 

320 

LGI1I -BAR 

c LSI  121)  *8** 

WRITE  I IS, *301  LC 
LGI1I-PLUS 
LSI  1211-PLUS 
IF  U8.50.0I  CO  TO  350 
18-0 

WRITS  I IS, 4*0)  LO 


N 1>1 
N 152 
N 153 
N 15* 
N 155 
N 156 
N 15T 
N 158 
N 159 
N 1*0 
N 1*1 
N 1*2 
N 1*3 
N 1** 
N 1*5 
N 166 
N 1*7 
N 1*8 
N 1*9 
N 170 
N 171 
N 172 
N 173 
N 17* 
N 175 
N 17* 
N 177 
N 178 
N 179 
N 180 
N 181 
N 182 
N 183 
N 18* 
N 185 
N 186 
N 187 
N 188 
N 189 
N 190 
N 191 
N 192 
N 193 
N 19* 
N 195 
N 196 
N 197 
N 198 
N 199 
N 200 
N 201 
N 202 
N 203 
N 20* 
N 205 
N 20* 
N 207 
N 208 
N 209 
N 210 
N 211 
N 212 
N 213 
N 21* 
N 215 
N 21* 
N 217 
N 218 
N 219 
N 220 
N 221 
M 222 
N 223 
N 22* 
N 225 
N 22* 
N 227 
N 228 
N 229 
N 230 


j 


Uk 


330 


€ 


I 


C 


€ 


340 

350 

360 

C 

C 

370 

3SO 


390 


W 

403 

4X0 

420 


43  0 
440 
450 
460 


00  340  1-1*121  N 

LDI  1 1 » )L  ANK  N 

CONTINUE  N 

CONTINUE  N 

N 

BOTTDM  UNDERLINE  N 

01  370  1-1,121  N 

LGI  II-  9 INUS  N 

01  380  1—1*121 ,12  N 

LGIl)-»LUS  N 

C ILL  SH3RT7  IVL.FR.TX)  N 

•WIT®  (IS, 4501  YL.LG  N 

WRITE  I IS, 4601  XI  H 

continue  n 

CALL  WS  CO, NX, NY)  N 

RETURN  N 

N 

FORMAT  I///21X,3A3,60X,A3,  1X.G12.3, //23H  FUNCT I ON-CONSTANT- ,G12  N 
1.4///I  N 

F JR  NAT  I///21X.3A8, 60X, A 8, 1X,G12.3,//40H  BLANK  PAGE  FOR  CONTOR.  N 
1 XB,XT,YL,VU-,4G12.3///>  N 

FORMAT  in  , 7X  ,4M-*  ,3Ad,  7X,60A1,10X,A8,1X,F7.3//14X,8HNAXINUM-  N 

1,  IP,  01 2. 1,  3X«  8HN1 NINUN-,  3 12.  J,  3X  , 13HCONTOUR  ZER0-,G12.  3 , 3X  ,1  8HC0NT  N 
2 OUR  INCREMENT-, G12. 3, 3X, AB/7  U,3P,  UIF8.1,4X))  N 

FORMAT  I 8X , 1 21 A1 ) N 

FORMAT  ( IH-.7X.121A1I  N 

FORMAT  ( IX ,F  7.  1,12141)  N 

FORMAT  ( IX, lit  F8  .1.4X1//// ) N 

ENO 


231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 


c 


c 


< 


< 


c 


c 


SU8R1JTINE  COMB IN  ( NL I M E , I OEST , I CM ,NFR ON, I FROM, I CHAR  I 0 

DIMENSI3N  IFROMI NFROMI,  ICHA4(  HFADM)  0 

INTEGER  PLUS, BLANK, ONE  0 

OATA  PL  iS/lH*/ .BLANK/1H  / 0 

DIMENSION  LGI140),  LBI14JI,  IBEGI10I,  IEN0I10),  IFILEI10I  0 

OATA  D/1H-.139-1H  / 0 

C 0 

C COMBINES  NLlNE  PRINT  LINES  FROM  THE  BEGINNING  OF  STREAMS  I FROM ( NFROM 1 , 0 
C ICHARINF4H)  CHARACTERS  FROM  EACH.  REWINDS  STREAMS  AT  THE  END.  0 

c deals  correctly  with  overprinting  ♦ at  line  beginning,  other  beginning  o 

C CHARACTERS  ARE  TREATED  AS  BLACKS  EXCEPT  ON  IFROMtl).  0 

C OUTPUTS  ICH  CHARACTERS  PER  LINE,  TO  STREAM  IDEST.  0 

C NFROM  AND  ICIAR  REDUCED  IF  NECESSARY  TO  GIVE  LESS  THAN  ICH  CHARACTERS.  0 
C IBM  STREAMS  //GO.FTXXFOOl  DO  UNI  T«  SY  SO  A,  SP  ACE  • I T RK , 1 1 , 5)  ) , 0 

C //  DC8-(RECFM-fba,BLKSI?£-1330,LRECL-133I  0 

C 0 

C 0 

C INITIALISE.  0 

N-O  0 

DO  3D  J-l, NFROM  0 

I3EGIJI-N-1  0 

N-N+ICHARIJ)  0 

IF  I (.LT.ICHI  GO  TO  10  0 

NFROM  » J 0 

N-  ICH  0 

ICHARI JI-ICH-IBEGIJ )♦!  0 

10  IENDIJI-N  0 

IFILEIJI-l  0 

K-l FROM ( J)  0 

IB-I3EGIJ)  0 

IE-IENCIJI  0 

REWIND  < 0 

C 0 

READ  14,20)  LM,ILGI1),1-IB,1E)  0 

IF  IJ.cl.l)  LL-LM  0 

C 0 

20  FORMAT  I 140A1)  0 

30  CONTINUE  0 

IC-N  0 

C 0 

C MAIN  LOOP  0 

00  93  ILINE-l.NL INE  0 

WRITE  I I OEST  ,20)  LL  , (LG  ( I ) ,1 -1 , 1 C)  0 

LL-8L A (<  0 

00  83  J-l, NFROM  0 

IF  (IFILE(J).EQ.OI  GO  TD  60  0 

IE-  IEMOI J ) 0 

I B-t  35G( J)  0 

K-IFRINCJ)  0 


10 


20 

30 


c 


•♦5 


( 


40 


c 


c 


READ  K, 20, END-601  LN,  I LSI  1 1 , I- I B,  I E I 
IF  ILH. IE. PLUS)  GO  TO  50 

WRITE  ( I OEST  (20)  (LB(I).I"1,H).  (LGI  I > , I • I 8, IE  I 
GO  TO  40 
50  IF  IJ.iJ.ll  LL-LM 

G'J  T 3 5) 

40  03  TO  !•  16  . IE 

TO  LGIII-3LANK 

IFILEI J)»0 

•o  continue 

90  CDNTINJ* 

C ENO  HA IN  LOO*. 

C 

00  100  J» l ,MFROH 
K-IFRONIJ) 

REWINO  < 

100  CONTINUE 

RETURN 
ENO 


0 50 

0 51 

0 52 

0 53 

0 54 

0 55 

0 56 

0 57 

0 58 

0 59 

0 60 
0 61 
0 62 
0 63 

0 64 

0 65 

0 66 
0 67 


I 

I 


c 


c 


SUBROUTINE  80  (A.II.JJI 

C •••  APPLY  STH NETRV  BOUND  AR Y C0N0IT13NS  FOR  ARRAY  PLOTTING  **• 
OMENS  UN  Al  I I.JJI 
III"! l-l 
JJl-JJ-1 
00  10  W.JJl 
AI1,UI"A<2,JI 
10  Al II, JI>AIII1, J) 

00  20  I « 1 . 1 ( 

All ,1 1 " A I 1,21 
20  Al I.J J)«AI  I, JJ  II 

RETURN 
ENO 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


c 


€ 


C 


( 


c 


c 


SUBROUTINE  YPLOT  (A,N,M,XY,X5,XT,CH, AL,AR,NVTK,NVCPT,NHTK,NHCPT, 
1 IS  I 
C 

C •••••  JlYN'S  YPLOT  PLOTTER  ROUTINE  *••••• 

C PLOTS  N FUNCTIONS  OF  X,  OOWN  T NE  PRINTER  PAGE. 

C X INCREASES  U«  THE  PAGE  FROH  Xi  TO  XT. 

C AIJ.KI  IS  THE  VALUE  OF  THE  K»TH  FUNCTION  AT  XVIU),  U«1  TO  N. 

C XV I J I is  VONJTONIC  INCREASING. 

C THE  LEFT  NARSIN  REPRESENTS  THE  VALUE  ALIK). 

C THE  RIGHT  RAISIN  REPRESENTS  THE  VALUE  ARIX). 

C THE  K*TH  FUNCTION  IS  PLOTTED  USING  THE  CHARACTER  CHIO. 

C THE  PLOTTING  AREA  HAS  NVTK  TICX  NARKS  VERTICALLY,  AN9 
C NVCPT  CHARACTERS  PER  TICK. 

c the  plotting  area  has  nhtk  tick  narks  horizontally,  and 
C NMCPT  CHARACTERS  PER  TICK. 

C THE  TOTAL  CHARACTERS  HORIZONTALLY  IS  NHT <"NHCPT , AN"  IS  LESS  THAN  132. 

c 

C BLANK  C UN non  is  USED  FOR  WORKING  SPACE  I 266  WORDS)  AND  SO  PASSED 
c ARRAYS  SHIULJ  NOT  BE  IN  BLANK  COHRJN. 

C 

IRPLICIT  REAL *81 A-H.O-Z I 
CUMMIN  LGI  13  3,  2) 

REAL  A 

0 I HENS  ION  AIN, Ml,  XV  IN)  , CHIN),  AL  M ) , ARM) 

DATA  9LANK/1 H /, MI NUS/ 1 H-Y ,*LUS/ 1H"/ , BAR/ 1 Hi / 

INTEGER  CM, BLANK, PLUS, BAR 

|H"NHTK»NHCPT 

I V"NVTK*NVCPT 

|H2"|H*2 

IVI"IV»t 

N4"N-1 

C WRITE  IIS.IDI 

ILP-2 

LGI  1,11  "BLANK 
L 'll!. 21-PLUS 
00  120  I"I,IV1 

IF  II.N-.l.ANO.I.NE.IVU  33  TO  30 
OJ  10  J"2, IH2 
LGI J, 1 1 "MINUS 
10  LGI J, 2 1 "BLANK 

00  23  U"2,IH2, NMCPT 
2D  LGI  J,  II  "PLUS 

60  TO  U 

>0  33  40  K-1.1L* 


0 

Q 

0 

0 

0 

0 

w 

0 

0 

0 

0 

Q 

0 

0 

0 

0 

Q 

0 

0 

0 

0 

Q 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1 

2 

3 

4 

5 

6 
7 
9 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 
IN 
20 
21 
22 

23 

24 

25 

26 

27 

28 
2° 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 


It 


A 


c 


f 


c 


c 


c 


c 


3?  VO  J-2.IH2 
43  LSI J,K)«8LAN< 

LSI IH2, ll-BAR 
L',(  2,  ll-BAR 
I HH-I 2-IH2 1/2 

IF  Il/2*2.E3.I.AN0.NHT</2*2.E.1.NHTK1  LSI IHH, 1 ) -8 AR 
53  COST  IN  JE 

ILP-1 

X-XT-I 1-11*1 XT-XBI/IV 
00  60  J-l.N1 

IF  (IX-XVI  Jl  l-IX-XVI  J-lll.LE.DI  SO  TO  TO 
60  CONTJNJE 

IF  ((X-XV(2)I/|X-XV(1>>.'ST.1|  J-2 
TO  J-MAX0IJ.2I 

J-NINOI J.N-2I 
IX- J-2 

0}  100  J-l,M 

c 

E-0 

00  90  <-1.4 

D-AIIX-K.JI 

OO  80  L-1.4 

IF  IL.EQ.Kl  SO  TO  80 

D-O-IX-XVI  IX-LII/IXVIIX-KI-XVI  IX  -L  ) I 

bo  continue 

90  E-E-0 

K-2 .5- ( c-AL  C J I l/IARI JI-ALI  Jl  l-IH 
IF  (K.jT.lH2.0R.K.LT.2)  SO  TO  100 
IL-1 

IF  UGlK.ll.NE.BLANKI  11-2 
LSI  X.  ID-CHI  J) 

UP-1AX0IILP.1U 
100  CONTINUE 

00  120  L-l.ILP 

IF  (I I-tl/NVCPT-NVCPT.EQ.I-l.AND.L.EQ.ll  00  TO  110 
WAITE  (IS. 1301  ILGIK.LI  ,K-l,H2) 

CO  TO  120 

110  LSI IH2 . ll-PLUS 

LSI  2.  ll-PLUS 

WAITE  (IS. 1401  X, ILGIK.LI, K-2, IH2I 
120  CONTINUE 

130  FORMAT  IA1.7X, 133A1 I 

140  FORNAT  ( IX ,F 7 . 1 , 133 A1 I 

RETURN 
ENO 


0 45 

Q 46 
0 47 

Q 48 
0 49 

0 50 

0 51 

0 52 

0 53 

0 54 

0 55 

Q 56 
Q 57 
0 58 

0 59 

0 60 
Q 61 
0 62 
0 63 

0 64 

Q 65 
0 66 
Q 67 
0 68 
0 69 

0 70 

0 71 

0 72 

0 73 

0 74 

0 75 

Q 76 
0 77 

0 78 

Q T9 
Q 80 
0 81 
0 82 
0 83 

0 84 

0 85 

0 86 
0 87 
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€ 


< 


C 


< 


SUBROUTINE  XPlOT  IF, XV. NX. I SEP , NL I N , XL . XU, TI TLE , SUBT  IT.VARII  R 

C R 

C •••  GLYN'S  XPLOT  PLOTTER  ROJTINE  **•*•  R 

C Fin  IS  VALUE  AT  XVI  tl.  MITH  XVI  1 1 NONOTONIC  INCREASING.  R 

C THE  DIMENSIONS  ARE  NX.  R 

C OOMAIN  XL  TO  XU  SPLIT  INTO  ISEP  PASES  HORIZONTALLY  FOR  HIGH  RESOLUTION  R 
C APPROX! NATELT  NLIN  LINES  ARE  OAAWN  VERTICALLY  ON  THE  PRINTER  PAGE.  R 

C TITLE  IS  24  CIARACTERS  AT  8 PER  WORD.  R 

C SUIT  IT  IS  3 CHARACTERS,  E.G.  -TINE  ■ • , WITH  VAR|  THE  TIME  VALUE.  R 

; MAIN  PROGRAM  SHOULD  CONTAIN  — R 

C COMMON  /TXTC/I TXTI 601, DAY  R 

C REA0I5,  II  ITXT  R 

Cl  FORNAT ( 6 0 A1 1 R 

C CALL  IOAXOAY)  R 

C — TO  READ  FJRTHER  TEXT  OFF  A CARO  AND  SET  DATE  IN  DAY.  R 

C BLANK  COMMON  IS  USED  FOR  WORKING  SPACE  (1568  WORDSI  AND  SO  PASSED  R 

C ARRAYS  SHOULD  NOT  BE  IN  BLANK  COMMON.  R 

C R 

IMPLICIT  REAL-81 A-H.0-2 I R 

INTE3E1  Jil 3,1211, JTI3, 1211  R 

REAL  G.8  R 

DIMENSION  FI  NX  I • GI721),  I HI  12 1 > • XI I ID  * XVINXI  R 

I NT  EG  E I KK  (31  R 

EQUIVALENCE  IKKI1I.KK1I,  (KKI2I.KK2I,  (KKI3I.KK3I  R 

I ITE5ET  012,2.21  .EOU.APOS,  JNO,  3L  ANK.M!  NUS,  PLUS,  B AR  ,00T  ,LG(  121 1 R 
D MENS  ION  TITLEI  31  R 

DIMENSION  ISTRMI10I  R 

COMMON  /TXTC/  ITXT 1 601 • OAY  R 

COMMON  J.JB.JT.IH  R 

DATA  Q/ IHM  , IN9 , 1FH,  1 HJ  » 1 HP  » 1 HI  ,1HL,1H  /.E0U/1H-/  R 

DATA  AP IS /1H* / , UNO/ 1H_/ ,J0T/1M,/,8AR/1H|/  R 

DATA  9LANK/1 M /, MI NUS/1H-V .PLUS/ 1N-V  R 

DATA  ISTRN/6, 21, 31, 4, 11,41, 42,43,22, 23/, NSTRM/10V  R 


1 

2 

3 

4 

5 

6 

7 

a 

9 
0 
.1 
2 
3 
.4 

5 

6 
.7 

8 
.9 
!0 
tl 
!2 
!3 
!4 
!S 
16 
!7 
8 
!9 

10 
>1 
12 
13 


liT 


< 


13 

23 

33 

C 

C 


40 


so 


c 

60 


73 

• 0 


C 


40 

130 

110 


NST-1 

R 

34 

IF  INLIN.EQ.-1  1 NST-NSTR1 

R 

35 

IF  INLIN.EQ.-ll  NLIN-16 

R 

36 

03  3)0  IS-l.NST 

R 

37 

IS*tSTUINSI 

R 

38 

XG-IXU-XLI/ISEP 

R 

39 

03  330  1SP-1.ISEP 

R 

40 

XB-XL4X3RISP-XG 

P 

41 

00  1)  1*1,121 

R 

42 

ton » *ii  nos 

R 

43 

03  23  1*1.121, 12 

R 

44 

lgi  1 1 **lus 

R 

45 

03  33  1*1.11 

R 

46 

XI I I > ■ < 44X  G*  1 1 — 1 1/  10 

R 

47 

R 

48 

NN»NX- I 

R 

49 

F I >40  14XI1U1,  NINIPRJH,  4N0  C ONT  O JR  ORIGIN  XNO  1NCRENENT 

R 

50 

C XLL  11X1IN  IF .XV.NX.FHX.F <N,F4d,X1X,X1N,X4BI 

R 

51 

IF  1 1 . £4*1 F1X-FNNI .GT»04dSIF1XI*343S( FHN ) I GO  TO  40 

R 

52 

WRITE  CIS. 3401  TITLE. SU6TI T.V4RI .FNX 

R 

53 

RETURN 

R 

54 

CO XT  IN Ji 

P 

55 

4OLJG10IIFNX-FNNI  /NLlNI  ♦.  1S3S 

R 

56 

1*4 

R 

57 

IF  I4.LT.0I  l-l-l 

R 

58 

B*I 

R 

59 

C*4-I 

R 

60 

F INC* l ).  **8 

R 

61 

IF  (C.LT..1S33  1 U TO  50 

R 

62 

FINC-IJ.  5 

R 

63 

IF  1C. LT.. 30131  GO  TO  SO 

R 

64 

F INC* 1 ) . **8*  2 

R 

65 

IF  IC.LT.. 4)161  GO  TO  S3 

R 

66 

F INC*  1 3 .••B*3 

R 

67 

IF  IC.LT..6S3SI  GO  TO  53 

R 

68 

F 1 NC*  1 J.  **B*  5 

R 

69 

IF  1 C .LT • • 84  90  I GO  TO  S3 

R 

TO 

F INC-1  3.  •*B*T.  5 

R 

71 

CONTINUE 

R 

72 

FINC«FINC*4 

n 

73 

l-FH/F  INC-.  09375* 1030 

R 

74 

1*1-1303 

R 

75 

for»i*finc 

R 

76 

X-IF1X-F0RI/F1 NC*1 . 15625 

R 

77 

N*N*4-  1 

R 

78 

FINC-FlNC/4 
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